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Abstract: As the Higgs boson properties settle, the constraints on the Standard Model 
extensions tighten. We consider all possible new fermions that can couple to the Higgs, 
inspecting sets of up to four chiral multiplets. We confront them with direct collider 
searches, electroweak precision tests, and current knowledge of the Higgs couplings. The 
focus is on scenarios that may depart from the decoupling limit of very large masses and 
vanishing mixing, as they offer the best prospects for detection. We identify exotic chiral 
families that may receive a mass from the Higgs only, still in agreement with the /177 
signal strength. A mixing 6 between the Standard Model and non-chiral fermions induces 
order 0 2 deviations in the Higgs couplings. The mixing can be as large as 6 ~ 0.5 in 
case of custodial protection of the Z couplings or accidental cancellation in the oblique 
parameters. We also notice some intriguing effects for much smaller values of 6, especially 
in the lepton sector. Our survey includes a number of unconventional pairs of vector-like 
and Majorana fermions coupled through the Higgs, that may induce order one corrections 
to the Higgs radiative couplings. We single out the regions of parameters where /177 and 
hgg are unaffected, while the / 17 Z signal strength is significantly modified, turning a few 
times larger than in the Standard Model in two cases. The second run of the LHC will 
effectively test most of these scenarios. 
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1 Introduction 


When the Large Hadron Collider (LHC) began its data taking, possible extensions of the 
Standard Model (SM) at the TeV scale were already severely constrained: electroweak 
(EW) precision measurements accurately confirmed the structure of the gauge sector [1—3], 
a number of flavour violating observables showed no significant deviation from the SM 
predictions [4, 5], all direct searches of non-standard particles at LEP and Tevatron gave 
null results [2, 6, 7]. After the first run of the LHC, the lower bounds on the masses of new 
particles increased substantially [8, 9]. The crucial discovery of the Higgs boson [10, 11] 
and the measurement of some of its properties [12-14] supported the minimal realisation 
of electroweak symmetry breaking (EWSB), as predicted by the SM. Thus, the room for 
SM extensions further compressed. 

In this phase, it is essential to reassess the possibilities still open for non-standard 
physics close to the TeV scale. In this paper, we focus on new spin-1/2 degrees of freedom. 
In particular, we will assume that the scalar and gauge sector is the SM one, with one 
standard Higgs doublet. In general, additional dynamics in the EWSB sector may well be 
present, including corrections to the Higgs boson couplings to the electroweak gauge bosons, 
as well as additional scalar or vector states: here we do not consider these possibilities, 
assuming that such dynamics takes place at sufficiently higher scale, or it is sufficiently 
weakly coupled to the SM. In other words, we will study effective field theories containing 
the SM degrees of freedom plus additional fermions only, being agnostic on the ultraviolet 
completion at higher energy. 

When wondering what fundamental fermions exist in Nature, one may notice that the 
fermion held content of the SM appears whimsical in some respects. Each fermion family is 
a set of five chiral multiplets, II, e#, qL, ur and cLr, whose gauge quantum numbers are not 
explained within the SM, with the remarkable property to be anomaly-free. The number of 
families, three, is unexplained too. All SM fermions are massless before EWSB and, when 
the Higgs develops a vacuum expectation value (vev), they acquire a mass proportional 
to their Yukawa coupling to the Higgs doublet. The structure of the Yukawa coupling 
matrices is not determined by the SM symmetries either. Of course, some of these issues 
may find a convincing interpretation in very high energy extensions of the SM, such as 
grand unification or flavour theories, but in this paper we will take a phenomenological 
point of view and centre on the TeV scale only. 

Additional chiral fermions that are massless before EWSB are definitely worth to look 
for, as their mass is bound to the TeV scale; the classical example of a chiral fourth family, 
with the same held content as the SM ones, was ruled out long ago [15], but we will show 
that more exotic possibilities exist. On the other hand, chiral fermions that transform in 
a real representation or form vector-like pairs, with respect to the SM gauge group, can 
acquire a mass before EWSB. While in general such mass can be much larger than the EW 
scale, in a number of well-motivated extensions of the SM there are new, fermionic degrees 
of freedom close to the TeV scale. Here are some familiar examples: 

• Non-zero neutrino masses may be generated through the mixing with heavier fermions, 
typically sterile neutrinos. The latter may have a mass close to the EW scale. How- 
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ever very different scenarios are possible, spanning from the eV scale to the grand 
unification scale. 

• The dark matter energy density may be carried by new, weakly interacting fermions. 
If thermally produced, their mass should be close to the EW scale. 

• If the quantum stability of the EW scale is guaranteed by supersymmetry broken 
at the TeV scale, the SM gauge bosons shall be accompanied by gauginos, and the 
scalar bosons by higgsinos. 

• If the weak scale is stabilised by dimensional transmutation, via a new strongly- 
coupled sector that condenses at the TeV scale, a number of composite spin-1/2 
resonances may be present in the low energy spectrum. For example, in the scenario 
of partial compositeness [16, 17], the SM chiral fermions, or at least the heaviest ones, 
are accompanied by vector-like composite partners with the same gauge quantum 
numbers. 

Given the diversity of phenomenological and theoretical motivations, and the wide- 
ranging discovery reach of the LHC, in section 2 we undergo a classification of the fermionic 
extensions of the SM, as general as possible. Theoretical consistency requires the absence of 
gauge anomalies. In addition, all new fermions should acquire a large mass to comply with 
null direct searches. As we want to explore the new constraints that materialised after the 
Higgs discovery, we limit ourselves to those fermions that interact with the Higgs doublet via 
Yukawa couplings. We provide the full list of SM extensions with these properties, formed 
by up to four new chiral multiplets, and comment on larger sets of new fermions. We note 
in passing that there may be alternative phenomenological motivations to study fermions 
that do not interact with the Higgs, for example to avoid the flavour problem altogether, 
or to look for generic dark matter candidates; some complementary classifications along 
this line can be found e.g. in Refs. [18-20]. 

The rest of the paper is dedicated to the phenomenology of the fermionic SM exten¬ 
sions, in particular to identify the regions of parameters that survive to three broad classes 
of constraints: (i) EW precision tests; (ii) collider direct searches; (iii) Higgs boson cou¬ 
plings. Our purpose is to provide a comprehensive, comparative description of all possible 
sets of fermions, that are presently allowed and may have an observable effect at the second 
run of the LHC. Such analysis has several limitations that one should keep in mind: 

• We compute only the leading order corrections to the Higgs and gauge boson couplings 
due to the extra fermions, and we roughly extract the collider bounds on their masses 
and couplings from the available experimental papers. A precision analysis would 
require a dedicated study for each given set of new fermions (and it is already available 
for several specific cases). 

• We assume that the new fermions do not mix with the first and second SM families, 
in order to avoid the strong constraints coming from flavour observables (tree-level 
flavour changing neutral currents are absent). Indeed, the mixing with the third 
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family is sufficient to characterise the corrections to the EW and Higgs observables, 
that are our main subject of interest. In addition, the Higgs couplings to the light 
families are presently unconstrained. Note that the mixing with the third family can 
still induce flavour violating processes at one loop, especially in H-meson decays or 
oscillations, when the top or bottom quarks mix with new fermions. However, the 
corrections are suppressed with respect to the SM by the masses of the new fermions 
and their mixing with the SM ones. The constraints are typically mild, but in some 
cases may be complementary to those discussed in this paper, see e.g. Refs. [21-23]. 

• The new fermions are not supposed to form an ultraviolet complete theory, consistent 
up to a scale much larger than the TeV. Therefore, we do not impose constraints 
coming from the coupling evolution at high energies, such as vacuum stability, absence 
of Landau poles, or gauge unification. We have no pretension to determine the full 
theory. 

• We do not restrict the possible SM extensions using cosmological considerations, 
that rely in most cases on specific assumptions on the early Universe evolution. For 
example, we do not impose bounds on the relic abundance of the new fermions, based 
on the assumption of an initial thermal abundance. 

As a matter of fact, these points can be addressed only in a model-dependent manner. In 
specific, well-motivated scenarios, it would be worth to perform more precise computations 
and include additional constraints from the other sectors of the theory. 

In section 3 we discuss purely chiral sets of fermions, that is, fermions that are massless 
before EWSB. Fermions with an EW-invariant mass, that is, either a Majorana or a vector¬ 
like mass term, are discussed in section 4, if they are colourless, and in section 5, if they 
are coloured. With a little abuse of terminology, we will call ‘leptons’ all the colourless 
fermions, even when they do not mix with the SM leptons, and ‘quarks’ all the coloured 
ones, even when they do not transform in the fundamental representation of the colour 
group SU( 3) c . Finally, in section 6 we recapitulate the most interesting results of our 
analysis. 

For each sets of new fermions, we were confronted with the need to compute EW 
precision observables and Higgs couplings. Thus, we took the opportunity to collect all 
the relevant formulas in the appendices, that generalise well-known results to the case of 
a generic fermionic extension of the SM. In appendix A, we present the fermion-gauge 
boson couplings, the corrections to the S and T parameters, as well as to the Zff vertex. 
In appendix B, we discuss the fermion-Higgs boson couplings, both the tree-level and the 
loop-induced ones, and we briefly summarise the present experimental constraints on the 
Higgs couplings. 

2 Minimal fermionic extensions of the SM 

Let us consider the extension of the SM by additional fermions, classified according to their 
transformation under the SM gauge group SU(3) C x SU(2) w x C/(l)y, whose irreducible 
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representations can be denoted by (R c , R W ,Y). The SM extension is defined by the most 
general renormalizable Lagrangian involving the SM fields and a given set of extra chiral 
fermion multiplets. 

We wish each set of new fermions (i) to be phenomenologically viable, (ii) to be theo¬ 
retically self-consistent, and (iii) to modify the Higgs couplings. This leads to the following 
series of requirements: 

(i) No massless fermions after EWSB, except for the three SM neutrinos and gauge 
singlets. Indeed, massless fermions are phenomenologically forbidden, unless they 
have no colour ( R c = 1), no electric charge (Q = T 3 + Y = 0), and no coupling 
to the Z-boson (T 3 — tan 2 6 W Y = 0). The latter two conditions imply Y = T 3 = 0. 
These conditions would allow for a massless neutral component in the chiral multiplet 
(1, R w , 0) with R w odd, however the gauge symmetries permit a Majorana mass term 
for such a multiplet. 

(ii) No SM gauge anomalies. The fermionic extensions of the SM under consideration are 
intended as effective theories valid up to the multi-TeV scale, therefore they should 
cancel all SM anomalies self-consistently. (Extra fermions much heavier than the EW 
scale play no role in the anomaly cancellation, since they form vector-like pairs with 
respect to the SM gauge group.) Since the SM field content is anomaly-free by itself, 
the anomaly-cancellation conditions must be imposed on the set of new fermions only. 

As we require the absence of massless coloured states, the new fermions form a 
(reducible) real representation of SU( 3) c , therefore the 517(3) c -cubic anomaly is au¬ 
tomatically vanishing. Denoting the new fermion representations by (R C i, R W i,Yi), 
for i = 1 ,..., n, the remaining anomaly-cancellation conditions read 

SU(3) C - SU(3) C - U(l) Y : ££=1 N wi C{R ci )Y = 0 , 

SU{2) w - SU(2) W - U(1) Y : £” =1 N ci C{R wi )Y = 0 , 

U(1) Y - U(1) Y ~ U (l)y : £” =1 N ci N wi Y? = 0 , 1 j 

grav - gray - U (l)y : Y!i=\ N ci N wi Yi = 0 , 

where N w = dim(i?,,„), N c = dim(i? c ) - this notation is redundant for 517(2) but not 
for 51/(3) - and the index C(R) of a given representation is defined by Tr(TjjTjj) = 
C(R)5 ab , with the index of the fundamental conventionally normalised to C(N ) = 1/2 
for SU(N). In the case of 517(2) one has C(R W ) = N U1 (Ay, — 1)/12. In the case of 
517(3) each representation R c is characterised by two integer Dynkin labels ( 01 , 02 ) 
with a* > 0, and one has N c = (1 + oi)(l + a 2 )(l + oi/2 + 02 / 2 ) and C(R c ) = 
N c {a\ + 3ai + 01 O 2 + 3o 2 + 02)/24. 

Additionally, the SU(2) w gauge group has a global anomaly, that cancels only when 
the sum £” =1 N c iC(R w i ) is an integer number [24], Note that C(R W ) is half-integer 
for N w = 2 + 4n, n = 0 , 1 , 2 ,..., and integer in all other cases. As for the previous 
anomalies, this condition must be satisfied by the fermions below the multi-TeV 
scale. (Heavier fermions, decoupled from the EW scale, necessarily give an integer 
contribution to the sum: only an even number of multiplets with N w even can acquire 
a vector-like mass.) 
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(iii) Non-zero corrections to the Higgs boson couplings. This corresponds to consider only 
new fermions with a Yukawa coupling to the SM Higgs doublet. More precisely, any 
subset of new fermions that satisfies the requirements (i) and (ii) by itself - any subset 
with no massless states nor anomalies - should have a non-zero Yukawa coupling to 
the Higgs. Otherwise, such subset would interact with the SM only through gauge 
interactions; strictly speaking, it can still affect the Higgs boson couplings at the 
two-loop level, but here we neglect such small effects. 

We stress that the three requirements above are independent, in the sense that none is 
automatically implied by the others. In particular: vector-like fermions are automati¬ 
cally massive and anomaly-free, but they may not couple to the Higgs doublet; chiral 
fermions that have non-zero masses, such as an extra family of quarks, can be anomalous; 
an anomaly-free set of fermions, such as zero-hypercharge fermions, may contain some 
massless components. 

In the following we classify the sets of n chiral fermions that satisfy the requirements 
(i), (ii) and (iii), for n = 1, 2, 3,4, and we briefly comment on larger sets. For convenience, 
we will mark with the symbol □ each viable set that is identified. 

2.1 One multiplet 

If we add to the SM only one new chiral fermion ~ (R c , R W ,Y), the only possibility to 
avoid massless states with non-zero SM gauge charges is the presence of a Majorana mass 
term that requires R c = R c , Y = 0 and N w odd. Such multiplet is anomaly-free. 

The additional requirement to couple to the Higgs doublet, H ~ (1,2,1/2), restricts the 
possibilities to R c = 1 and R w = 1 or 3, that is, 

□ N ~ (1,1,0) or E~ (1,3,0). (2.2) 

In both cases a Yukawa coupling is allowed among the new fermion, H and the SM lepton 
doublet l ~ (1, 2, —1/2). 

Since N (X) forms by itself a self-consistent extension of the SM that modifies the 
Higgs couplings, n replicas of N (and/or X) also define a set of new fermions satisfying 
all our criteria. We will analyse their phenomenology in section 4.1. Of course, there may 
also be consistent sets of n new fermions that are partly formed by replicas of IV or E, and 
partly by different multiplets, as we will see in the following sections. 

2.2 Two multiplets 

Let us classify the possible pairs of chiral fermions ipi and that can be added consistently 
to the SM and that modify the Higgs couplings. The fermion t/i can satisfy all requirements 
without -02 only if it transforms as (1,1, 0) or (1, 3,0), as shown in section 2.1. In this case 
the three obvious possibilities are 

□ N\ + N 2 , Ei + E 2 , IV + X . (2-3) 

For all other representations, a coupling between ipi and ^2 is necessary: either there is 
a vector-like mass term or a Yukawa coupling {H). The latter possibility 
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leads to an inconsistent mass spectrum: one has N w i = N w 2 + 1 (or vice versa) and the 
unbalanced component of ipi has T 3 7 ^ 0, therefore it cannot be massless. Then, either 
ipi admits a Majorana mass term or it couples to a SM fermion too. One can check that, 
in both cases, either another unwanted massless state is left, or a SM neutrino acquires a 
large mass too. An alternative way to exclude the case of the Yukawa coupling 
is to solve the anomaly system ( 2 . 1 ) for n = 2 : one obtains N c 1 = N c 2 , C(R c 1 ) = C(R c 2 ), 
R w ± = R w 2 and Y\ = —Y 2 . Thus, we conclude that the two chiral fermions should form a 
vector-like pair, 


□ ipi ~(R c ,Rw,Y), ^2 ~ (R c , Rw, —Y) ■ (2.4) 


In order to modify the Higgs couplings, at least one among tpi and ^2 should have a 
Yukawa coupling with a SM fermion and the Higgs doublet. We take conventionally all 
chiral fermions to be left-handed. A SM family is formed by l ~ (1,2,—1/2), e c ~ (1,1,1), 
q ~ (3, 2,1/6), u c ~ (3,1, —2/3), and d c ~ (3,1,1/3). In order to have a Yukawa coupling 
with these representations, the new fermions should transform under SU( 3) c either as 
singlets or triplets. The former mix with SM leptons and can be called vector-like leptons 
(VLLs), the latter mix with SM quarks, hence the name vector-like quarks (VLQs). Under 
SU(2) w they can transform as singlets, doublets or triplets. All possible vector-like fermions 
with a Yukawa coupling to the SM fermions are listed in Table 1. 

To analyse the phenomenology of vector-like fermion multiplets, it is useful to name 
the components with different electric charge Q. The possible components of the VLLs 
have charges Q(N ) = 0, Q(E) = — 1 and Q(F) = —2. Then, the self-conjugate leptons N 
and E and the four VLLs can be written as 


N, 


E = 


( E °\ 

N ; E, 

w 



A 



(2.5) 


After EWSB, N, E and E c can mix with the SM leptons v, e and e c , respectively, while 
F does not mix with the SM. We will discuss the phenomenology of these VLLs in section 
4.2. 

The possible components of the VLQs have charges Q{X) = 5/3, Q(T ) = 2/3, Q(B) = 
— 1/3 and Q(Y) = —4/3. They are embedded in seven possible VLQ multiplets, 



( 2 . 6 ) 

After EWSB, T, T c , B and B c can mix with the SM quarks t, t c , b and b c , respectively. 
On the contrary, the components X, Y and their conjugate do not mix with the SM. We 
will discuss the phenomenology of these VLQs in section 5.1. 
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IpLl 

1pL2 

coupling to 

E~ (1,1,—1) 

E c ~ (1,1,+1) 

l,e c 


L c ~{ 1 ,2,+^) 

l,e c 

A ~ ( 1 , 2 ,—|) 

A c ~ (1, 2, +|) 

e c 

(1,3,-1) 

A c ~ (1, 3, +1) 

l 

T~ (3,l,+§) 

T c ~(3,l,-§) 

q,u c 

(3,1,-|) 

H c ~(3,l,+i) 

q,d c 

X T ~ (3,2,+|) 

Xj. ~ (3,2, — g) 

u c 

Q~(3,2,+|) 

Q c ~( 3,2,-|) 

q, d c , u c 

Yb ~ (3,2, — |) 

YJ~(3,2,+|) 

d c 

Xq ~ (3,3, + |) 

Xq ~ (3, 3, — |) 

q 

y q ~ (3,3,-|) 

Y4~(3,3,+i) 

q 


Table 1. Vector-like pairs of left-handed chiral fermions, that provide a consistent extension of the 
SM and modify the Higgs boson couplings. 

2.3 Three multiplets 

Let us classify the possible sets of three chiral fermions ^1,2,3 that can be added to the SM 
consistently with the requirements of section 2 . Of course, there is the trivial possibility 
to combine smaller sets that are already consistent on their own: 

□ Three copies of N ~ (1,1,0) and/or of S ~ (1, 3,0). 

□ A vector-like fermion from Table 1 plus one copy of N or X. As the latter couples 
to the SM lepton doublet l , there may be a non-trivial interplay with the VLLs E, L 
and A, that also couple to l. We will discuss this case in section 4.4. 

To explore all the other possibilities, note that there are two patterns for the colour repre¬ 
sentations of V’1,2,3) that guarantee the absence of massless coloured states: 

• R c \ = R cl plus a vector-like pair R c 2 = R c 3 7 ^ R, c \ . The only choice allowing for 
couplings to the Higgs is R c 1 = 1 and R c 2 = 3. In this case the two subsets ipi and 
ip 2,3 do not interact with each other, therefore each should be a consistent extension 
of the SM by itself, reducing to the trivial possibilities already listed above. 

• Rci = R C 2 = Red- These can be either real or complex representations. One may 
have considered three real representations of SU(3) C not all equal to each other; in 
this case, however, it is not possible to allow Yukawa couplings to the Higgs and to 
give a mass to all coloured components, at the same time. 
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In the only non-trivial case, R c \ = R c2 = R c 3 or permutations, the anomaly conditions in 
Eq. (2.1) reduce to 

' ^ 1+^*2 + ^3 = 0 , 

< N wl y\ + N w2 Y 2 + N w3 Y 3 = 0 , (2.7) 

k N w iYf + N w2 Yi + N w3 Y 3 3 = 0 . 

The possible solutions of this system (with N w i positive integers) are 

(a) Yi = 0 and N W i arbitrary, for i = 1, 2, 3. 

(b) Y\ = 0 with N w 1 arbitrary, Y 2 = —Y 3 0 with N w2 = N w 3 , or permutations. 


(c) For any pair of integer numbers n, m > 0, 


N w i — 1 T n , N w 2 
Y\ 7 ^ 0 arbitrary, Y 2 = 


2 + n + m , N w3 = 3 + 2 n + m , 

y V JV»1 +2 ~ Kl 

1 JV„ 2 (JV2,3 - JV2 l2 ) ’ 3 1 JV^s JVj 3 - ivj 2 ’ 


or permutations. 


( 2 . 8 ) 


Let us consider these solutions in turn, to analyse whether they can satisfy the other 
requirements of section 2 . 

(a) In order to couple to the Higgs boson, one needs R c i = 1 and R w i = 1 or 3 for 
z = 1,2,3, that is three copies of N or E: a trivial possibility already considered. 


(b) Suppose first that il> 2 and ip 3 form a vector-like pair, that is, R c \ = R c2 = R c 3 . 
But, to avoid unpaired complex representations of SU( 3) c , one needs i? c i to be real. 
Barring the trivial cases where 'i/q and ^ 2,3 form self-consistent extensions of the SM 
separately, the necessary condition to modify the Higgs couplings is to allow for a 
Yukawa coupling between -i/’i and ^ 2 , 3 - This leads to 

□ Vh ~ {RciR W i 0 ), ^2 ~ (R c ,R w ±i,-), 1P3 ~ (R c ,Rw± 1, — -), R C = R C - 

(2.9) 

The choice of R w is arbitrary up to the SU(2) w global anomaly: for N w = 2 + An, 
one needs N c to be even. The colourless case R c = 1 will be discussed in section 4.4. 
The coloured case R c = 8, 27,... will be discussed in section 5.3. 

Next, suppose that if) 2 and 1/23 do not form a vector-like pair, that is, a complex 
representation R c \ = R c2 = R c3 . In order to have the same number of colour- 
conjugate representations one needs N w 1 = 2N w2 . A Yukawa coupling is also needed 
among 'i/q and ^ 2,3 to provide masses, so the only possibility is 

□ ipi ~ (R c , 2 , 0 ) , ijj 2 ~ (R c , 1 , -) , iji 3 ~ (R c , 1 , --) , R C ^R C - ( 2 - 10 ) 

The SU(2) w global anomaly further requires that N c must be even. Note that this 
is the minimal, consistent set of chiral fermions that has no vector-like mass terms, 
rather it acquires a mass from the Yukawa couplings only. The phenomenology of 
purely chiral fermions is discussed in section 3. 
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(c) For (n, m) / (0, 0), one has N w j, > 4. Then, -^3 does not couple to the SM nor to -01, 
and even the possible Yukawa coupling to if ;2 cannot provide a mass to all the N w s 
components of ^3. Therefore, let us take (n, m) = (0,0), that implies 

Vh ~ {Rc, > ^2 ~ (-R c , 2, --Y) , ^3 ~ (R c ,3,-Y) . (2.11) 

By requiring an equal number of components with opposite electric charge, one finds 
Y = ±1/2 or Y = ±1/6. In both cases one can check that some components of 
the new fermions remain massless, therefore no consistent SM extension of this type 
exists. 

2.4 Four multiplets 

In the previous sections we derived the list of all the consistent sets of n new chiral fermions, 
with n < 3, discussing in detail how to implement the requirements of section 2. Here we 
provide the complete list for n = 4, without displaying the lengthy and involved analysis 
needed to prove this result. 

First of all, there are a number of possibilities to combine smaller subsets of new 
fermions that are already consistent by themselves. It is worth to list them for bookkeeping 
and to point out those combinations with special phenomenological relevance: 

□ Four copies of N and/or £. 

□ Two copies of N and/or £ plus a vector-like fermion from Table 1. In particular, 
the set (IV, £, L, L c ) corresponds to the neutralinos and charginos of the minimal 
supersymmetric SM: bino, wino and the two higgsinos. This case is discussed in 
section 4.4. 

□ Two vector-like fermions Ti and T2 from Table 1. A non-trivial interplay occurs 
when Ti and ^2 couple both to a given SM fermion, as indicated in the last column 
of Table 1, and/or when there is a Yukawa coupling between Ti and T2: this happens 
for Bor A with L or A; T or Xq with Xt or Q, B or Yq with Q or Yb- In models of 
partial compositeness, a SM fermion acquires its mass by mixing with two vector-like 
composite fermions, with the same quantum numbers as the SM left- and right- 
handed components: Q and T for the top quark, Q and B for the bottom quark, L 
and E for the tau lepton. This case is discussed in section 4.3 for leptons and 5.2 for 
quarks. 

□ One copy of N or £ and a set of three fermions from Eq. (2.9) or Eq. (2.10). A 
non-trivial interplay occurs in the case (1,3,0) + (1,5,0) + (1,4,1/2) + (1,4, —1/2), 
discussed in section 4.4. 

Let us come to the consistent sets of four fermions that are not the union of two smaller 
self-consistent sets. We found that non-trivial solutions are possible only when the four 
colour representations R c i are all equal or conjugate to each other. After all requirements 
of section 2 are taken into account, only two possible patterns emerge: 
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• For arbitrary R c and R w , a viable set of four multiplets is provided by 

□ ( R c , Rw — 1, 0 ) + (R c , Rw + 1, 0 ) + (R c , Rw, 1/2) + (Rc, Rw, ~1 /2) , (2-12) 

with one exception: if N w is odd, then either C(R W — 1) or C(R W ± 1) is half-integer, 
therefore one needs N c even to cancel the global SU(2) w anomaly. 

The case R c = 1 can be described as two Majorana leptons plus a vector-like lepton 
(see section 4.4). The case R c = R c 7 ! 1 is the analogue for coloured fermions (section 
5.3). Finally, the case R c 7 ! R c is purely chiral, with no masses before EWSB (section 
3). 

• For arbitrary R c , R w and Y, a viable set of four multiplets is 

□ (Rc, Rw, Y) + (R c , R w , —F ) + (Rc, Rw + 1, F — 1/2) + (Rc, Rw + 1, ~Y + 1 / 2 ) . 

(2.13) 

For the first time we encounter a pattern where the hypercharges of the new fermions 
are not determined uniquely. 

The case R c = 1 corresponds to two VLLs (section 4.3), except when Y = 0 with N w 
odd, or Y = 1/2 with N w even: then, one has two Majorana leptons plus one VLL 
(section 4.4). The case R c 7 ! 1 corresponds to two VLQs (section 5.2), except when 
R c = R c and Y = 0 with N w odd, or Y = 1/2 with N w even: then, one has two 
Majorana quarks plus one VLQ (section 5.3). 

We found that all other sets of four multiplets relevant for Higgs couplings are not viable: 
either some component remains massless, or a gauge anomaly is present. 

2.5 Larger sets of new fermions 

We do not attempt a general classification for n > 5 new chiral fermions. On the one 
hand, the general principles and the different phenomenological possibilities are already 
well illustrated by more minimal sets of fermions. On the other hand, a detailed analysis is 
worth only in the context of a specific, well-motivated theory beyond the SM. Here we shall 
mention some prominent examples that have been extensively studied, to situate them in 
the context of our classification. 

• We have shown that there are two sets of purely chiral fermions, displayed in Eq. (2.10) 
and in Eq. (2.12), formed by three and four multiplets, respectively. The more tra¬ 
ditional chiral extension of the SM is a fourth family, formed by the five multiplets 
q ' L , t' R , b' R , l' L and t' r . It was already excluded at the time of LEP, because the Z 
invisible width forbids a fourth massless neutrino, but it could be rescued adding a 
sixth multiplet, a sterile neutrino v' R . It is by now excluded by the measurement of 
the Higgs boson couplings [25], as we will review at the end of section 3. 

• The minimal supersymmetric extensions of the SM predicts fermionic partners for the 
gauge bosons and for the two Higgs doublets. This amounts to five chiral multiplets: a 
bino ~ ( 1 , 1 , 0 ), a wino ~ ( 1 ,3, 0 ), two higgsinos ~ ( 1 , 2 , ± 1 / 2 ) and a gluino r\_/ ( 8 , 1 , 0 ). 
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The latter does not enter in our classification, since it does not couple to the Higgs 
doublet. Concerning the other four multiplets, supersymmetry restricts the possible 
couplings among them and to the SM, therefore it corresponds to a special case in the 


parameter space of the SM extension by the set (N, E, L, L c ). We will briefly discuss 


the related phenomenology in section 4.4. Of course, our purely fermionic extension 
of the SM corresponds to the limit where the scalar supersymmetric partners are 
significantly heavier than neutralinos and charginos. 

• Another scenario addressing the hierarchy problem is compositeness. An effective 
way to couple the SM fermions to a composite Higgs doublet amounts to a partial 
fermion compositeness: each SM chiral fermion mixes with a composite vector-like 
fermion with the same quantum numbers. Thus, to induce a Yukawa coupling among 
two SM fermions and the composite Higgs one needs two vector-like fermions. There¬ 
fore, a SM extension by four chiral multiplets is suitable to study this mechanism for 
one SM Yukawa coupling at a time. Indeed, partial compositeness corresponds to a 
special subspace of parameters, because the symmetries of composite models restrict 
the couplings of the new fermions and of the SM ones. We will briefly discuss the phe¬ 
nomenology of r-compositeness in section 4.4 and the case of b and f-compositeness 
in section 5.2. Of course, realistic models of partial compositeness require more than 
two vector-like fermions, e.g. to induce the Yukawa couplings of all the heavy, third 
family fermions. The interplay between the Higgs and composite vector-like fermions 
is studied in detail e.g. in Refs. [26-30]. 

3 Phenomenology of new chiral fermions 

In this section we study purely chiral sets of new fermions, that is to say, one cannot write 
any fermion mass term before EWSB, thus their masses are generated by the Higgs vev 
only. These sets do not contain vector-like pairs of chiral multiplets, that would admit a 
Dirac mass, nor multiplets in the representations (R C ,R W , 0) with R c = R c and R w odd, 
that would admit a Majorana mass. A purely chiral set, consistent with the requirements of 
section 2, constitutes a new fermion ‘family’, very much analogue to the three SM families. 

We identified two classes of purely chiral sets, formed by three and four multiplets 
respectively, displayed in Eq. (2.10) and Eq. (2.12). We will discuss the phenomenology of 
these two classes in some detail. In the last part of the section, we will investigate whether 
larger chiral sets of fermions may be compatible with present Higgs data. 

• Three chiral multiplets. The only consistent ‘family’ formed by three chiral multiplets 
is 



(3.1) 


the smallest viable representation being R c = 6. The Yukawa interactions are 


+ Aia^iL-fD/YR + h.c. . 


(3.2) 
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Here and in the rest of the paper we do not display the obvious kinetic terms, that must 
be added for each new fermion. After EWSB one is left with two mass eigenstates ±12 and 
F 13 in the same colour representation R Cl with charge Q = ±1/2 and mass m\2 = Ai2w/\/2 
and m \3 = Ai 3 u/\/ 2 , respectively. 

The lightest new fermion is stable and forms hadrons with exotic charges, that are 
constrained to be heavier than several hundreds of GeVs. Indeed, the searches for R- 
hadrons [31, 32] assume the existence of a stable stop or sbottom (scalar with R c = 3), or 
of a stable gluino (fermion with R c = 8). In the latter case one finds m Rc= s > 1320 GeV 
[33], and we expect a similar (stronger) bound for R c = 6 (larger colour representations), 
as these fermions are pair-produced through their coupling to gluons. To roughly estimate 
the limit on a stable sextet, we rescale the gluino bound by computing the ratio between 
the sextet and octet production cross section at tree level. Taking into account the different 
colour contractions and the interference of the s, t and u-channel, we obtain a lower bound 
m Rc= 6 > 1400 GeV. Further discussions and references on the constraints on R-hadrons 
can be found e.g. in Refs. [20, 34]. This limit from direct searches already leads to some 
tension with the perturbativity upper bound, 77112,13 <C {4n)v/y/2 — 2.2 TeV. 

The contribution of ± 12,13 to the oblique parameters reads 


N c 
S ~ — 
6vr 


T ~ 


N r 


167rs?„c?„m 


2 1 2 

m 12 + m 13 
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m 12 m 13 
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m \ 2 
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12 


— m 


13 


m 


13 


(3.3) 


Note that for (mi2 — 77113) —> 0 one finds T oc (mi2 — mis) 2 , because in the degenerate limit 
the custodial symmetry is restored: for A 12 = A 23 Eq. (3.2) has a global SU(2)r symmetry 
with ip 2 R and i/jsr transforming as a doublet. The value of S can lie within the 3cr ellipse, 
but only if 1V C = 6 and T ~ 0.3 at the same time (see figure A. 2). This can be achieved 
for (mi2 + mis)/2 ~ 1500 GeV and (777-12 — 77113) — 50 GeV. Thus, this ‘family’ of new 
fermions is marginally compatible with direct searches and EW precision tests. 

After the Higgs discovery, however, one can definitely exclude such a set of new 
fermions. As they are heavier than the Higgs boson and do not mix with the SM fermions, 
the Higgs decay width at tree-level are unchanged. However, a huge deviation occurs in 
the loop-induced Higgs coupling to gluons, 


R 


99 


°(gg -± h) 

&sM{gg -± h) 


~ [1 ± 4 C{R c )f . 


(3.4) 


where we have taken the limit , m\ 2 , 7nf 3 3> tti|/ 4 in the loop form factor (see appendix 
B.2). Even the smallest possible colour representation has C(6) = 5/2, leading to a huge 
Rgy = 121, totally incompatible with LHC data. 


• Four chiral multiplets. Let us move to the only consistent ‘family’ formed by four 
chiral multiplets, 

, i/j 4R ~ 

(3.5) 



i>lL ~ {Rc, Rw-l, 0), V±L ~ (Rc, Rw + l, 0 ), 7/3 R 


F c , Rw 7 2 
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with R c R c to prevent vector-like or Majorana mass terms, and with N c even if N w is 
odd, to prevent a global SU(2) w anomaly. The allowed Yukawa interactions are 

— Cy = + Xi^IlHi/j^r + \ 23 'i/J 2 lH'iJj 3 R + A24V’2L#V’4B + h.c. . (3.6) 

where the explicit SU(2) W contractions are defined in Eq. (B.4). After EWSB the compo¬ 
nents of l combine with those of 03 , 4 /j to forms 2 N w mass eigenstates: one with charge 
Q = N w /2\ two mixed states with Q = N w /2 — 1,..., — N w /2 + l] one with Q = —N w /2. As 
they have (half-)integer charges, they do not mix with the SM quarks and the lightest state 
is stable and hadronises, with collider bounds above 1 TeV, analogue to those discussed 
above. The discussion of oblique parameters is also similar to the previous case: one can 
be marginally consistent with data, by choosing the parameters to realise an approximate 
custodial protection. 

The way to definitely exclude this set of chiral fermions is, once again, their contribu¬ 
tion to the Higgs boson coupling to gluons. Note that each of the 2 N w mass eigenstates 
belongs to the same colour representation R c and must be (much) heavier than the top 
quark, therefore one finds 

R gg ~ [1 + 2(2N w )C(Rc)} 2 . (3.7) 

Even in the minimal case with N w = 2 and C( 3) = 1/2, one finds a very large R gg — 25, 
incompatible with the LHC Higgs data. 


• Larger sets of chiral multiplets. Let us ask the question whether we can exclude any 
set of purely chiral fermions, even when it is formed by more than four multiplets. Indeed, 
the Higgs coupling to gluons implies that any new chiral fermion should be colourless, 
because even the minimal set of chiral coloured fermions, formed by a weak doublet and 
two weak singlets with R c = 3, leads to a large R r f g n — 9. This is not compatible with 
the range currently allowed by global fits, 0.5 < R gg < 1.8 at 99 % C.L. [35] (see appendix 
B.3 for details). In particular, in this way one can exclude [25] a fourth SM family, formed 
by the six multiplets q ' T , t' R , b' R , l' L , t r and v' R . Recall that the sterile neutrino is required 
to avoid an additional massless neutrino, that is forbidden by the Z invisible width; then 
this set of fermions is not purely chiral, but one may invoke a lepton number symmetry to 
forbid the sterile neutrino Majorana mass. Let us remark that coloured chiral fermions may 
be allowed in the case of extended Higgs sectors, not considered in the present paper. For 
example, adding an Higgs triplet, it is possible to rescue the fourth family [36]. Another 
example is provided by coloured chiral fermions receiving their mass from a second Higgs 
doublet [37]. 

Coming to colourless chiral fermions, some of the mass eigenstates are necessarily 
charged and thus contribute to the Higgs width to photons as 


R -y ry - 


I 477 i All 
I^SAf ' sine// 

\All | 2 

I "Asm I 


477 ^ 

new — 




(3.8) 


where i ? 77 is defined in Eq. (B.41), the SM amplitude is A”p M — —6.5 and the sum runs 
over the new fermion mass eigenstates. Note that, to derive Eq. (3.8) from Eqs. (B.21) 
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and (B.22), we took (i) 2m* mh, that is accurate enough, even though the lower bounds 
on heavy charged lepton masses are weaker than those on coloured particles; (ii) Higgs- 
fermion couplings yt = rrii/v and fji = 0 , that is the case for purely chiral fermions, because 
their mass matrices are proportional to the Higgs vev v, see Eq. (B. 6 ). Note also that the 
result is independent from potential mixing between the new fermions and the SM leptons. 
The presently allowed range is 0.5 < < 1.9 at 99 % C.L. [35]. For n chiral multiplets 

( R w i,Yj ,), one finds 


477 ~ 2 SAST'r,2 = 2 ST' + 12Y i 2 ~ !) 

' /x new — 2 / , / , h/fc 2 / , -^2 


(3.9) 


i= 1 k= 1 i =1 

We added an overall factor 1/2 to take into account that each massive fermion is formed by 
two chiral components. Equivalently, one may take the sum only over fermions of a given 
chirality. For example, a minimal set is formed by a weak doublet 2^ of hypercharge Y 
paired with two weak singlets (1 + 1)#, giving a contribution Anew = 2(1 + 4Y 2 )/3 > 2/3. 
The next-to-minimal set (2 + 2)l paired with (1 + 1 + 1 + l)/j implies AnJw > 4/3, that 
is still allowed by the present constraint on i? 77 , while for example (3 ± 2)l gives already 
AZlw > 10/3, that is almost excluded. Since the SM amplitude has opposite sign w.r.t. 
the one of new fermions, one can also envisage the contrived possibility of a large AZew ~ 

_9 477 ^ 1 Q 

^SM — 10 ‘ 

This shows that there are purely chiral sets of n fermions that satisfy the 77 con¬ 
straint. However, we have also shown before that no set exists for n < 4, that satisfies 
the consistency requirements of section 2. It is non-trivial to check whether purely chiral 
sets with n > 4 could be consistent. Consider for example the case of two weak doublets 
plus four weak singlets. In order for all components to receive a mass from the Yukawa 
couplings to the Higgs doublet, the hypercharges must be chosen as 


(2, hi) , (2,Y 2 ) , (l,-Hi + 0 , (l,-Hi-£) , (l ,-H 2 + £) , (l,-Y 2 -0, 

(3.10) 

in the convention where all multiplets have the same chirality. The absence of anomalies 
requires Eq. (2.1) to hold, and this leads to Y\ = —Y 2 = Y. As a consequence, three 
vector-like mass terms are allowed and such set of fermions does not qualify as purely 
chiral. 

At this point one should recall that, in this paper, we took the point of view that all 
mass terms and interactions allowed by the gauge symmetries are present. In alternative, 
one can easily introduce some global symmetry to forbid possible vector-like mass terms, 
thus imposing by hand that a given set of fermions is chiral. For example in Eq. (3.10) take 
a U{ 1) symmetry with charge +1 for doublets and —1 for singlets. If one takes this point 
of view, the set of fermions in Eq. (3.10) qualifies as the minimal still viable set of purely 
chiral fermions. Indeed, besides being consistent with all the requirements of section 2, it 
can be compatible with direct searches, EW precision tests, and constraints from the Higgs 
couplings. 

Concerning direct collider searches, as the new leptons have charges Q = Y ±1/2, they 
do not mix with the SM leptons (except for |Y| = 1/2 or 3/2) and the lightest state is 
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\Q\ 

1/3 

2/3 

1 

2 

3 

4 

5 

6 

7 

8 

bound in GeVs 

200 

480 

574 

685 

752 

793 

796 

781 

757 

715 


Table 2. The 95 % C.L. lower bounds on the mass of heavy stable leptons, from the CMS 
collaboration [32]. The production is assumed to occur through the Drell-Yan process only. Limits 
are obtained for SU(2) w singlets, but they remain similar in general [20]. The ATLAS collaboration 
obtains comparable but less stringent limits in the range 2 < \Q\ < 6 [38]. For larger values of Q 
see Ref. [39]. 

stable. There are severe bounds on the number density of such charged relics [20], but they 
depend on cosmological assumptions: e.g. for a reheating temperature below their mass, 
they were never produced in the early Universe. Limits on heavy stable leptons are of 
the order of a few hundreds of GeVs, and are displayed in Table 2 for some representative 
values of Q. For a review on heavy stable particles see Ref. [34], 

The contributions to the S and T parameters of a fermion ‘family’ formed by one 
doublet and two singlets are given in Eqs. (A.24) and (A.25). In the present case we have 
two such ‘families’ with N c = 1 and opposite hypercharges ±T, and one can easily lie 
within the 3<r ellipse of Fig. 10. For example in the custodial limit where the two mass 
eigenstates of each ‘family’ are degenerate, one finds T ~ 0 and S — 0.1. Concerning the 
Higgs couplings, one finds AZlw = 4(1 + 4T 2 )/3, that lies in the allowed range of R 77 for 
|y| < 0.3 and 1.4 < \Y\ < 1.6. 

Finally, it is interesting to compare Eq. (3.9) with the analogue amplitude for the Higgs 
boson coupling to 'yZ. For n chiral multiplets ( R w i,Yi)i one finds 

jpZ ^ 2 g T^k - s^Qk _ 2 y-y - 12T - 2 tan 2 9 W - 1) (3 11) 

^ i= 1 k =1 ° w ^ i= 1 

where we used Eq. (B.31) particularised to the case of chiral fermions, in the same way we 
did above for the 77 case. Note that, when summing over the mass eigenstates of equal 
charge Q, the mixing matrices disappear from the Z couplings in Eq. (A. 6 ), therefore one 
reduces to a sum over the interaction eigenstates. For example, the set in Eq. (3.10) gives 
Al Z ew ~ 2 [1 - (1 + 8 T 2 ) tan 2 9 W \ /3. 

4 Phenomenology of non-chiral leptons 

In this section we discuss new colourless fermions, which admit either a Majorana or a 
vector-like mass term before EWSB. 

4.1 Majorana leptons 

Let us consider leptons that admit a Majorana mass term. The latter requires a vanishing 
hypercharge, Y = 0. The two possibilities relevant for the Higgs couplings are sterile neu¬ 
trinos N ~ (1,1,0), and weak triplets E (1,3,0). 
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• Sterile neutrinos. In the case of one sterile neutrino, the SM Lagrangian is extended 
to £ = Csm + £tv, where 

Cn = -lLa^NaHN R - -N^M n N r + h.c. . (4.1) 

Only one linear combination of the three active neutrinos v Ra couples to Nr, and 
since we are not concerned with flavour issues we will drop the index a in the following. 
After EWSB, v R and Np mix and combine into two Majorana fermions and Uf, with 
definite masses, m Ul < m Uh ] active neutrino mass searches imply m Ul < 1 eV. If one takes 
the limit m Vt -C rrih, the type I seesaw mechanism is realised: m Ul — \ 2 N v 2 /{2 Mjv) and 
m Uh — Mjy. The Higgs boson couplings to vm and VhVh are proportional to m Vl /v , and the 
coupling to viVh is proportional to y/m Ul m Uh /v, leading to a negligibly small decay width, 
T(/i —> vivh) < m 2 h /{&TTV 2 )m Ul < 10 -8 MeV. Also, one can check that the decay widths of 
the Z-boson receive negligible corrections, always proportional to m u Jv. For m Uh above 
the EW scale, one can compute the contribution to the S and T parameters, that turns 
out to be suppressed by the tiny ratio m Ul / m Uh . 

In the case of two or more sterile neutrinos IVj, in most of the parameter space the 
arguments above apply to each separately: either the Majorana mass Mi is as small as 
the eV scale, or the active-sterile mixing 6i = Ajv* (y/\/2)/M^. is suppressed, ~ 

m vi < 1 eV. In either case the corrections to the Higgs and Z/W -boson couplings are 
tiny. The only exception occurs when much larger &i are tuned among each other, in order 
for the Ni contributions to the light neutrino mass to cancel. Consider for simplicity two 
sterile neutrinos N± At leading order in the mixing angles 9i one has 

m Vl ~ 1 0 2 M Nl + 0 2 M N2 1 < 1 eV . (4.2) 

The two summands have a physical relative phase, therefore they can be orders of magni¬ 
tude larger than rn Ul , if there is a strong cancellation between the two: the active-sterile 
mixing can be large, no matter how large the sterile masses m Vhl 2 are. Even though this 
scenario requires a severe tuning of parameters to lead to observable effects, it may be justi¬ 
fied by some symmetry. For example, in the so-called inverse seesaw model [40-42] (see also 
Ref. [43]), the lepton number symmetry U(1)l is broken by a small mass parameter, and 
the cancellation occurs naturally in the limit where this parameter goes to zero. Therefore, 
it is worth to analyse the phenomenological consequences of a large active-sterile mixing: 
both Higgs couplings and EW gauge boson couplings can be significantly modified. 

Consider first the neutrino mass eigenstates vi, z^i, in the regime m Vl <C rn Vhl 2 <C 
mh, mz■ One finds that the decay widths of the Higgs boson can be significantly modified, 
in particular T(h -t uiVhi) — m h ml h .\9i\ 2 /(8irv 2 ) and T(h -» v hi v hi ) ~ m h ml hi \0i\ A /(Anv 2 ). 
These rates can be easily as large as the total SM Higgs width, T? M ~ 4 MeV, therefore 
the LHC experiments already constrain 0* and m Vhi . Note that both invisible and visible 
decay channels are affected, since Vhi ,2 decay not only into light neutrinos, but also into SM 
particles e.g. via virtual W-bosons. Detailed analyses of the parameter space and of various 
constraints can be found e.g. in Refs. [44-48]. Note that the Z-boson invisible width T" w , 
that is measured at the few per mil level, is not significantly affected for m Ufil 2 < 1 MeV, 
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with Vhi ,2 decaying mostly invisibly into three ui\ even in the presence of large mixing, only 
the active components of vi.h\,h 2 couple to the Z-boson, and one recovers the SM value 
of T 1 ™ once the sum over all neutrino pairs is taken. On the contrary, for larger m Uhl 2 
the heavy neutrinos mediate visible Z-decays, therefore T 1 ™ is depleted and a significant 
upper bound applies on \9i\. 

Consider next the complementary regime mh,mz < m u hl2 - this case the Higgs 
and Z decay width are not modified, but a significant active neutrino fraction in 1 /^ 2 can 
still have observable consequences. Direct searches of EW scale sterile neutrinos through 
their mixing with active neutrinos have been performed e.g. by ATLAS [49, 50] and CMS 
[51, 52], Here we would like to point out that the EW precision parameters S and T can also 
receive important corrections, that constrain the masses and mixing of the sterile neutrinos. 
To understand this quite surprising fact, that is generally overlooked, it is convenient to 
write the 3x3 neutrino mass matrix in the basis [ vl ,, N^, CVJ) 2 ) as 

M u = U*diag(m, Vl ,m Vhl ,m Vh2 )U t , (4.3) 


with U unitary and subject to the constraint 0 = {M v )n — Ui 2 m v hl + U*%m Vh2 , where 
we neglected the tiny m Vl . Then, the active neutrino fractions Uu contained in the mass 
eigenstates can be parametrized in full generality as follows: C /13 = 9 taken to be real, 
C /12 = i6/y/rh with = rnu hl /m Uh2 and [C/n| 2 = 1 — 9 2 (l + rh)/rh- The 3cr lower bound on 
Tf v implies 6 2 {1 +rh)/rh < 0.015. We computed T and S using the formulas in appendix 
A, as a function of 6 , and rz = neglecting the mass of the SM leptons and 

including a symmetry factor 1/2 for loops of Majorana fermions. Here we report the result 
in some physically interesting limits: 


S = 1 /( 4 + 3 lo ® 7/f 


29 2 


( b ) rriu h2 » m„ hl = m z : T ~ 



87 Ts 2 cf„ m% 


3 - log 



g ^ 92 m "H2 

47r mz 


(4.4) 


In case (a), taking the maximal allowed value 9 2 iax — 0.007, the correction to T grows 
quadratically with the sterile neutrino mass: requiring to remain in the 3<r ellipse in the 
S — T plane (see Fig. 10), one finds the upper bound m Vh2 < 8.5 TeV. This sensitivity 
to very large scales is due to the significant fraction of the active neutrino in the heavy 
states; note that this non-decoupling effect requires a strong tuning among the two sterile 
neutrino parameters. In case ( 6 ), the active fraction in the heaviest sterile neutrino is 
rather 9 2 nax ~ 0.015(m^/'m^ 2 ), therefore T grows only logarithmically with m Uh2 , while S 
remains constant: one remains in the ellipse for m Uh2 as large as the Planck scale. 


• Weak triplets with zero hypercharge. In the case of a weak triplet Y,r (1,3,0), 
the SM Lagrangian is extended by 


= - 


- -Tr (Ej/A/eEr) + h.c. , 


(4.5) 
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where we adopted the matrix notation 


5 ^ = y/2 Z a R r a = ± 




^ - i^ R 


&R + &k ~& R 


( J_vO 
= [ ^ R 

V 


-s 


R 


_LyO 

V2^R, 


(4.6) 


and we normalised the triplet Yukawa coupling according to appendix B.l. After EWSB, 
the neutral component YP r and a linear combination of active neutrinos combine into two 
Majorana fermions v and So, in complete analogy to the case of v\ and ur discussed above. 
As usual, we will consider only the mixing with the third lepton family, taking An e ,S/i = 0 
and Asr = As- Indeed, flavour changing neutral current processes such as n —> ey or 
H —> 3e are strongly constrained [53, 54]. In the limit X^v <C Ms one realises the so-called 
type III seesaw mechanism: m u ~ A|u 2 /(6 Mv;) and ?ns 0 — Ms- The charged components 
(S^) c and E^ mix with the SM charged leptons tr and tr respectively, to form the mass 
eigenstates r and E _ . Then, three real parameters - the Yukawa coupling As, the mass 
Ms and the SM tau Yukawa coupling A r - determine the mass of four physical states, m v 
and M s o in the neutral sector, m T and M s - in the charged sector. The mixing angles 
for neutrinos and left-handed charged leptons are suppressed by yAn^/Afvo; the mixing of 
right-handed charged leptons receives an additional suppression by m T /Ms-. 

The couplings of So to the Z and Higgs bosons are exactly the same as the couplings of 
Vh discussed above. In particular, for vanishing m u all Z and Higgs couplings reduce to their 
SM values, therefore the corrections are negligibly small. At tree level, M s - — M s o also 
vanishes with , however it is well-known that at one loop weak interactions induce a mass 
split, M s - — M s o ~ 170 MeV (see e.g. Ref. [18]). Neglecting the tiny mixing angles, the 
only heavy lepton couplings are ZE + E _ and 1E + E“E°; other mixing-suppressed couplings 
are relevant for E-decays [55]. The contribution of E to the EW precision parameters S 
and T is vanishingly small, as EWSB is felt only through the mixing angles and through 
the loop-induced mass splitting among the E-components, and both are very small. 

Coming to direct searches, LEP looked for new charged leptons pair produced and 
decaying to Wu, setting a lower bound M s - > 100 GeV [56]. At the LHC heavy leptons 
are mostly pair-produced via Z* / 7 * — > E + E _ and W±* —> E^E 0 . The fraction of E that 
decays into each lepton flavour, b a = 6 a /(6 e + 0^ + 0 T ), characterises the final state. CMS 
[57] considered either b e = = b T = 1/3, b e = 1 or = 1, obtaining constraints in the 

range Ms > 180 — 210 GeV. The most stringent constraint comes from ATLAS [58], with 
Ms > 325 GeV for b e = 1 and Ms > 400 GeV for 5 /t = 1. We expect a weaker bound in 
the case b T = 1, that we assumed above. Other decay channels relevant for E searches at 
the LHC are discussed in Ref. [55], including displaced vertexes, as E becomes long-living 
in the limit of very small mixing. 

In the case of two or more lepton triplets E,;, the phenomenology is similar, except 
when the mixing between the SM leptons and the new leptons is not suppressed. As in 
the case of sterile neutrinos, this is possible only by severely tuning the Yukawa couplings 
of the various S ? ; to keep m u small. In the case of two triplets, the neutrino mass matrix 
is diagonalised as in Eq. (4.3), while the charge lepton mass matrix can be written as 
M e = ULdiag(m T , M^, M-£ 2 )U r . Neglecting m u and m T , the left-hand mixing matrix Ur 
coincides with the neutrino mixing matrix up to a y/2 factor: (Ul) 13 — \/2#, (Ur) 12 — 
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iy/29/y/rh and \(Ul)ii \ 2 — 1 - 29 2 (1 + r^/r^, with 77 = M Si /Me 2 , while the mixing 
angles in Ur are further suppressed by m T /M When 0 is large, the corrections to S and 
T may become significant as already discussed for sterile neutrinos. In addition, the new 
charged leptons, that are necessarily above the EW scale, could contribute significantly to 
h —> 77 , 7 Z. Before computing these corrections, one should notice that a strong upper 
bound on the mixing comes from the Z coupling to t + t~ . The LEP measurement of 
r(Z —> t + t~) [3] implies 29 2 (l + r^/r/j < 0.004 at 3o\ 

Let us describe in some detail the main corrections to R 77 and R^z, defined by 
Eq. (B.41). Similar analytic approximations could be used for the models analysed in 
the next sections as well. For the diphoton channel, using the results of Appendix B.2.2 
one finds 


R. 


-Arm + (\(U L )u\ - 1)A 1/2 (t t ) + \(U L )u\ A 1/2 (r Sl ) + |(C/l)i 3 | A 1 /2 (te 2 ) 


77 — 
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SM I 


Asm + 202 h ' A i/2{ T T, k 

• h 


1/177 I 

I-Asm I 


1 i0 2 1 + r h A l/ 2 (TEk) 

I'h 


All I ~ 

•^SM I 


> 0.998 , 


(4.7) 

where we took the maximal allowed values for the mixing angle and the form factor, 
A.i /2 (t s fc ) — 1.5 for — 100 GeV. Note that R 77 ~ // 77 , because the Higgs produc¬ 
tion rate and the total Higgs width are not significantly modified with respect to the 
SM. Thus the diphoton signal strength can be slightly reduced (the fermionic part of the 
amplitude slightly increases and interferes destructively with the IT-loops), but only at 
a few permil level. For the 7 Z channel the new physics contribution can be written as 
■Anew = .Ay- Z diag + A']f 0 ff-diag • Using the results of Appendix B.2.3, the loops involving a 
single mass eigenstate give 
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(4.8) 


where we took the Z couplings to the interaction eigenstates, thus neglecting corrections 
of higher order in the small mixing. The loops involving two mass eigenstates give 

\(Ul)ii\ 2 \{UL)ik\ 2 


A-Zfiff-diag 5 A= 2 , 
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/ = A 1 / 2 {tt, A r , rs fe , AeJ - z^^==B 1 / 2 (t t , A r ,rs fe , A Sfc ) 

V m T Mx k yJm T MY. k 


(4.9) 


Retaining only terms of order 9 2 and neglecting m T /Ms k , the rate relative to the SM can 
be written as 
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We neglected the B L / 2 term, as it interferes only with the very small imaginary part of 
the SM amplitude. The diagonal and off-diagonal form factors have comparable size, 
M/itokA-Zk) - !-3 and v / M Sfe /m T Ai/ 2 (T T , A r , rv fc , Asj.) — 1, where we took the large 
Ms k limit. Replacing the maximal allowed value for the mixing, we find R^z > 0.998, 
with a suppression at the few perrnil level, of the same order as for 


Finally, let us note that, with two or more sterile neutrinos Ni (or triplets E,) the CP 
symmetry can be broken. In general, the Higgs couplings to the fermion mass eigenstates 
are not real, and the off-diagonal couplings of the Z-boson can be complex as well. This 
does not modify any of the above results, because the CP violating effects vanish in the 
limit m Ul /m Uhk —> 0 (and m r /my; fe —> 0 in the triplet case), therefore they are subleading. 


4.2 One vector-like lepton 

Let us consider the addition to the SM of one vector-like lepton (VLL). The four different 
possibilities are a weak singlet E, a weak doublet L or A, a weak triplet A, whose charges 
are displayed in Table 1. As usual, we restrict ourself to mixing with the third SM family, 
i.e. with r and v T . The SM Lagrangian is extended by 


- £</, = \^l L Hip R + M^ip L ip R + h.c. , ip = E,A, (4.11) 

- £</> = A ^ l H{H)t r + M^ip L ip R + h.c. , ip = L( A) , (4.12) 

where the SU(2) w contractions are understood (see appendix B.l for details). In the case 
of E (L) one could write an additional mass term rieErtr (mi^LR), but such term can 
be removed by choosing conveniently the basis for the two fields Er and tr (II and Ll), 
that have identical charges. Thus, in each case there are only two real parameters: the 
vector-like mass M^ and the Yukawa coupling A^,; the mixing among the new leptons and 
the SM ones vanishes for vanishing A^,. There is no CP violation. 

The components of each multiplet ip are listed in Eq. (2.5). The doubly-charged 
component F does not mix as there is no SM counterpart with Q = 2, therefore m 2 F = M 2 . 
The <5 = 1 component E mixes with the SM r to form the two physical mass eigenstates 
t' and t. The mass matrix is given by 


M e 


( Xt V2 

VO ) 


ip = E,A , M e 


( 0 ) 

V A ;/; /,; J 


ip = L, A ■ (4.13) 


The SU(2) W Clebsch-Gordan coefficient is equal to one, except in the triplet case, 
= ^/l/3. The rotation to the mass basis can be parametrized as 


M e 


V L 



Vr 


Vl,r 


( C L,R Sl,R 
\—SL,R Cl,R 


(4.14) 


The triangular mass matrix structure of Eq. (4.13) implies some strict relations among the 
mixing angles and the mass eigenvalues. For the case of L and A, one finds 

771 777 777 / 

tan 9 l = —- tan Or < tan Or , s L = -t^sr , c L = c R . (4.15) 

m T f 
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For the case of E and A, the same relations hold with L v R. Note that direct searches 
of charged leptons at LEP [56] require m T > > 100 GeV , therefore one angle is at least 
two orders of magnitude smaller than the other. In the following we will refer only to the 
dominant mixing angle 9^ for each ip, dropping the subscript L,R on (co)sines. Note that 
m 2 , ~ A/|/c^ > A/ 2 ,. The neutral component N does not mix with the SM neutrino u T in 
the case of L. In the case of A there is mixing in the neutral sector, described by 


M v = 





(4.16) 


(4.17) 


with s 2 ~ 2s\/(l + s a) > s^. One neutrino remains massless, while the second acquires 
a mass m 2 , ~ (1 + s\)m 2 , > m 2 ,. In summary, the tree-level spectrum of heavy leptons 
satisfies 

Me < rry (ip = E) , mj? = A/a < r?y < Try (ip = A) , 
r?V = A II < 77v O’ = L) , mj? = A/a < m T > (ip = A) . 
with the mass splitting controlled by the mixing, Am 2 (ip) ~ s 2 A/ 2 . 

Let us briefly discuss the collider bounds on . In first approximation one can neglect 
the mass splitting. It is possible to recast some LHC multi-lepton searches to put bounds 
on VLLs. The limits on AZ^ strongly depend on the SM generation that couples to the 
heavy leptons. For couplings only to the third one and for the doublet L, Ref. [59] reports 
A/a > 280 GeV, while the LEP limit remains more constraining in the case of the singlet 
E, Me > 100 GeV. For the exotic doublet A with a doubly-charged component, Ref. [60] 
reports M\ > 320 GeV. To the best of our knowledge, no similar analysis is available for the 
triplet A. We expect a bound comparable or slightly stronger than to the one for A. These 
bounds only apply for promptly decaying particles. We will only consider this possibility, 
because heavy leptons become long-lived (cr > lm) for a tiny mixing — 10 -8 — 10 -9 , 
and the mixing suppresses all the deviations from the SM that we are interested in. More 
details on the collider phenomenology of A and A can be found in Refs. [61] and [62, 63], 
respectively. 

It is mandatory to require that the Yukawa coupling A^ lies in the perturbative regime, 
|A^| <S 4 - 7 T. This consistency requirement translates into an upper bound on the product 
of the heavy lepton mass and the mixing angle, 


I-VI 


V2 


77V 


Av -0 


V 


-S'ljj 


<C 47T 


(4.18) 


The perturbativity constraint on the SM Yukawa coupling A r is satisfied a fortiori. The 
non-zero couplings of the physical Higgs boson to the charged mass eigenstates, using the 
convention of Eq. (B.5), are given by 
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(4.19) 

where the plus (minus) sign holds in the case of E and A (L and A). In the case of A, 
there are also non-zero couplings to neutral leptons, 
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(4.20) 
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Important constraints come from the Z-decays into SM leptons. The couplings of the 
fermion mass eigenstates to the Z are defined in Eq. (A. 6 ). The couplings to the SM 
leptons tl, v t l and tr are modified if they mix with new leptons with a different weak 
isospin T 3 . Neglecting (m T /m^)-corrections, at tree level one finds 


R(Z 
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(gV T y + (g? T ) 


A \2 
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V,SM \ 
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A,SM\ 
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(4.21) 


where receive a correction of order s 2 with respect to the SM. The experimentally 
allowed range given in Eq. (A.28) implies an upper bound on the mixing angle, for any 
VLL: we find se,a 1$ 6.0 • 10 ~ 2 and sl,a 1$ 6.7 • 10~ 2 . In addition, in the case of A there is 
a correction to the Z-coupling to neutrinos, and thus to the Z-invisible width, 


i9uu) 2 + (9™) 


A \2 

lzizJ 
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6 6 Kdv'v Y + ygvv ) 


(4.22) 


The couplings gYl A receive a correction of order s 2 , that leads to a comparable limit sa 
8.2 • 10 -2 . Extracting the couplings of and B M from Eq. (A. 4), one can calculate 

the S and T parameters with the general formulas in appendix A.2. We find, at leading 
order in the mixing angle, 
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(4.23) 


where a^, and are numerical coefficients of order one. Taking into account the upper 
bound s^p < 0.06 from Z -A t + t~ , as well as the perturbativity bound from Eq. (4.18), 
(m T i/mz)sf < 10, we checked that S and T always lie in the allowed ellipse of Fig. 10. 

Coming to the Higgs boson signals, we first recall that all the dominant Higgs pro¬ 
duction channels at the LHC are not affected by the new leptons, as they leave the Higgs 
couplings to gluons and quarks unchanged. The total Higgs width also receives negligible 
corrections, as new leptons affect only the partial widths T(h -A a), for a = t + t~ , 77 , 7 ^. 
Therefore the Higgs signal is given by the ratio of partial widths in the model w.r.t. the 
SM, fj, a ~ R a . The tree-level Higgs decays are directly controlled by the couplings in 
Eq. (4.19) and Eq. (4.20), in particular 


R TT = (1 - 4 ) 2 > 0.99 , 


(4.24) 


where we used the bound from Z -A t + t~ . There is also the marginal possibility that 
the new leptons are lighter than the Higgs boson, thus opening the channels h -A tt' and 
h -A vv'. However, direct searches seem to allow the singlet E only to be sufficiently light. 
Using equation (B.7) and neglecting m T /m T i, we find 


T(h -A t + t' ) ~ 


c E t> E ,,l T' 
167TU 2 


m h 



< 0.2 MeV . 


(4.25) 


Note that both the couplings y TT ' and y TT > contribute equally to the decay width, see 
Eq. (4.19). We maximised the product s|.m 2 ; by taking se = 6 • 10 “ 2 and m T > = 100 
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GeV. As the SM total Higgs width is ~ 4.1 MeV, an enhancement of order 5% may be 
possible. Note that experimental searches at the LHC concentrated on h —> rr [64, 65], rp 
[ 66 , 67] and pp [ 68 , 69]. These channels are suppressed due to the small masses of the SM 
leptons, in contrast with the Higgs decays into a SM lepton plus a heavy lepton. It would 
be interesting to perform a dedicated search for this channel. 

For the photon-photon channel we find 
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(4.26) 


where the form factors are defined in appendix B.2.2. The addition of a VLL amounts to 
an additional r 7 -loop and a modified r-loop, that interfere destructively with the W-loops. 
Maximising the mixing and choosing m T i = 100 GeV (the form factor decreases for larger 
masses), we find ~ —1.9 • 10~ 3 , a perrnil reduction of the signal strength. For the 7 Z 
channel the relevant Z couplings, g^^ TT , T , T ,, receive corrections of order s ^ relatively to 
their unmixed values, as follows from Eq. (A. 6 ). At leading order in the small mixing and 
neglecting m T , we find 
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(4.27) 


where . is the isospin of the Q = — 1 component of the multiplet ip, and the plus 
(minus) sign in front of the off-diagonal term corresponds to the case ip = L (ip = E, A, A). 
As a consequence, the diagonal and off-diagonal terms always interfere destructively. The 
relative magnitude of the form factors is given below Eq. (4.10). The size of the correction 
change depending on the VLL under consideration, but it is always very small. The 
maximal deviation is obtained for A, with 5R^z — 1-3 • 10~ 3 . 


4.3 Two vector-like leptons (including r compositeness) 

Let us consider a SM extension by two VLLs. They may couple to each other by a Yukawa 
interaction or not. 


• Two VLLs not coupled to each other. In this case, each VLL must be a consistent 
extension of the SM by itself, therefore it should have the quantum numbers of E, L, 
A or A, that are displayed in Table 1. The six possible pairs of VLLs decoupled from 
each other are (ip, ip') = (E, E'), (L, L'), (A, A'), (A, A 7 ), (E, A) and(L, A). In the first four 
cases the additional mass term ipip' can be rotated away without loss of generality. The 
phenomenological effects are a trivial sum of those discussed in section 4.2 for a single VLL, 
with one noticeable exception. 

When L and A have the same Yukawa coupling to the SM and the same mass, the 
Lagrangian 

- -Cl,A = (Ll A l) f tr + [l l A L ^j + ^- c - (4.28) 
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preserves a global SU(2)l x SU(2)r symmetry. In this custodial limit the corrections to 
the T parameter vanish, and those to the coupling Ztrtr vanish as well [70]. This is the 
smallest set of tau custodians [71, 72]. A linear combination of the charge-one components 
of L and A, t" = (E^ — E^)/y/2, does not mix. The orthogonal combination, E = 
(E( l ) + E^)/y/ 2, mixes with the SM exactly as shown in Eqs. (4.13)-(4.15), to form the 
mass eigenstates t' and r. The spectrum reads m T n = rrt/v = mp = M^ < m T / ~ M^/cr. 
As discussed in section 4.2, direct searches already require all these states to be heavier 
than the Higgs boson. Thanks to the custodial symmetry, the Z couplings to leptons 
do not constrain the right-handed mixing sr between r and t'\ one has 5g// r = 0 and 
s 9tt = 4/ 2 = ( m r /M^) 2 Sr/ 2, that is negligibly small. The Zvv coupling is SM-like as 
well. The T parameter is almost SM-like as no additional sources of custodial breaking are 
introduced, and the correction to the S parameter is within the experimental range. 

The most stringent constraint on sr comes from R TT ~ (1 — s 2 r)‘ 2 ■ Using the 3cr lower 
bound R tt > 0.2 (see Table 4), one finds sr < 0.7. Indeed, the mixing can be large 
and reduce significantly the hrr coupling. As a consequence, the total Higgs width may 
be slightly reduced and, correspondingly, the signal strength for the other Higgs decay 
channels, defined in Eq. (B.40), may augment by a factor Tf t M /Th 1$ 1.04. The deviation 
in the 77 channel has the same form as in Eq. (4.26): imposing the constraint from h —> rr, 
one finds a lower bound /j 77 > 0.86, that is close to the present experimental sensitivity. 
Coming to the 7 Z channel, the loop involving both r and t' vanishes when one neglects 
the tiny left-handed mixing sl, because the custodial symmetry imposes g?', = 0. Then, 
Eq. (4.27) reduces to 

R lZ ^ 1 + 24tan 2 0 w Al,2{T ^ K,) , (4.29) 

\^sm\ 

with a maximal correction 5/x^z — 0.12. In the near future the increasing experimental 
precision on r tt can further constrain or eventually determine the mixing parameter sr in 
this custodial limit. 

• Two VLLs coupled to each other, not mixing with the SM fermions. Next, we 
have to discuss the case of two VLLs coupled through a Yukawa interaction. The most 
general assignment for their four chiral components is 

VhL, V’lii ~ (1) Rw, Y) , 1p2L, ^2R ~ (1, Rw + C Y + -) , (4.30) 

with the Lagrangian 

— + ^21^2lHi/>ir + + T/ 2 '02L'02/? + h.c. . (4.31) 

The four phases of A12, A21, M\ and M 2 cannot be all rotated away: one phase is physical 
and allows for CP violation. In the special case Y = 0 (Y + 1/2 = 0) and R w odd (even), 
one should add Majorana mass terms for iJjil,r (^2 l,r)- we postpone to section 4.4 the 
discussion of sets formed by one VLL plus Majorana leptons. For a few other values of 
R w and Y , displayed in Eq. (4.41), interaction terms between 1/I72 and the SM leptons 
are allowed and should be added to the Lagrangian. We discuss first the no-mixing case 
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and postpone to the end of this section the discussion of the mixing with the SM. In the 
absence of mixing, the lightest new lepton ipnght is stable, at least at the renormalizable 
level. If Y is an integer multiple of 1/2, ipu g ht may decay into a SM lepton through some 
higher dimensional operator. For all other values of Y, this state is absolutely stable and 
it has non-zero electric charge. Collider searches put a lower bound on the mass of stable 
heavy leptons as a function of their charge Q, see the discussion in section 3 and the limits 
in Table 2. The effect of the set of fermions in Eq. (4.30) on h —> 77 was studied e.g. in 
Ref. [73]. 

Let us begin by analysing the case R w = 1. There is one state with Q = Y + 1 and 
mass M 2 , and two states with Q = Y that mix, with mass matrix 
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l 777,21 M 2 
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(4.32) 


As My is the most general 2x2 matrix, it is useful to parametrize it in terms of the five 
physical parameters, 
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where 7771,2 are the real and positive masses of the eigenstates fi :2 , the mixing angles 6l 
and Or vary between 0 and 7r/2, and the CP violating phase ip varies between 0 and 27r. 
The only restriction comes from the perturbativity of the Yukawa couplings, 
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(4.34) 

These relations imply e.g. an upper bound on the masses for fixed values of the mixing 
angles. Vice versa, as the masses become larger and larger, the mixing angles vanish and 
the new fermions decouple from the EW scale. 

The Higgs boson couplings to fi t2 are directly obtained from Eq. (B.5). Taking for 
illustration the limit where fi and f 2 are mass-degenerate, one has |Mi| = IM 2 I = M, 
|mi 2 | = 1 77721 1 = H and mi = m 2 = M 2 + \j? = In this case the contribution to the 
amplitude for h —> 77 is 


•^h h = 1Q 2 F{6 l ,0rMAi /2 (77.) , F{0 l ,0rM = s 2 l c 2 r + c 2 l s 2 r - 2c L s L c R s R cos ip . 

(4.35) 

The CP-odd contribution vanishes because y\\ = —y 22 in the degenerate limit. The per¬ 
turbativity conditions in Eq. (4.34) reduce to F(0l,0r,tt) <C 87t 2 v 2 /m 2 ,. The interference 
with the SM is destructive as A"fjf 2 > 0. There are two allowed regions of parameters: 

(i) A SM-like region for small Q: the smallness of the charge ensures a small, negative 
departure from the SM. 

(ii) A cancellation region at large Q: for Q 2 ~ 4.8 /F(0l,0 r ,ip), the rate is accidentally 
close to the SM as A^ ~ —2AU rm . 
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Figure 1. Signal strengths ju 77 (red) and /r 7 z (blue) in the case of two VLLs ipi ~ (1,1, Q) and 
■02 ~ (1,2,(2 + 1/2), as a function of Q. We chose the following mass matrix parameters, defined 
by Eq. (4.33): mi = m 2 = 800 GeV, p = 0, and three values for the relevant mixing angle, 
9l — 9r = 7 t /8 (dotted), 7 r /10 (dashed), 7 t /12 (solid). The grey parts of the curves are excluded by 
S and T, see Eq. (4.36). The shaded horizontal band is the presently allowed range for at la 
(dark) and 3 a (light). 


This behaviour is illustrated in Fig. 1 for the CP conserving case p = 0, where F(0l, Or, 0) = 
sin 2 (0l — Or). Indeed the amplitude grows from ip = 0 to tp = n, as F{6l,0r,tt) = 
sin 2 (Or + Or). Further constraints on the parameters come from the EW precision tests. 
In the limit mz we find 
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(4.36) 

where M is the mass of the unmixed state with Q = Y + 1. Therefore, for fixed values of 
the mixing parameters, S and T constrain the charge Q, as illustrated in Fig. 1. 

It is interesting to analyse the value of /i 7 z in the allowed space of parameters. Using 
Eqs. (A. 6 ) and (B.31), the CP-even amplitude for h —> 7 Z is given by 


jO z 

^/i ,/2 


~2 ~F ( 0 l , Or, p) ( 4 l + 922) A1/2 ( Txf,, \ip) , 

Cm 


9ii F 922 2 ^Q s w ■ (4.37) 


Note that the loops involving both f\ and /2 vanish as they are proportional to the Higgs 
coupling 2 / 12 , or to the form factor Bi/ 2 defined in Eq. (B.39), and they both vanish for 
mi = m 2 - As the SM amplitude is negative, the new contribution interferes constructively 
as long as Q( 1 + 4s 2 , Q) > 0. The CP-odd amplitude f 2 is also non zero, because 
9ii T” 922i an( 4 ^ can b e sizeable for large values of sin tp. Let us distinguish the two regions 
of parameters allowed by /r 77 : 
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(i) In the SM-like region at small Q , varying (Qr — Or) we find —0.01 < Sp^z 1$ +0.08. 

(ii) In the fine-tuned region at large Q we find 2.5 < p^z 1$ 3.2 for tp = 0. This range 
slightly depends on the sign of Q , see Fig. 1. It is actually possible to obtain an even 
larger Pyz, while keeping // 77 close to one. For example, taking for simplicity sr = 0, 
the CP-odd amplitude reads 

A'jfj 2 =-4Q 2 c ls l tan 2 OwSiiupA^ij^,^) . (4.38) 

For sin ip of order one, this contribution becomes important and one can reach p^z — 
7. 


When one allows for the two mass eigenstates f\p to be non-degenerate, m\ < m 2 , the 
amplitudes for the diphoton channel become 
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where we made the approximation / 2 {t\ ) ~ ^ 1/2 (T 2 ), that is accurate for 4mf 2 
m^. For sufficiently large mass splitting the interference of AJJ ^ with the SM can be 
constructive. In the 7 Z channel the analytic form of the amplitude becomes more involved, 
in particular the loops involving both and fi are non-zero, and the interference with 
the SM strongly depends on the ratio mi/m 2 - One can tune the parameters to cancel the 
corrections to /r 77 in Eq. (4.39), e.g. taking <p = 0 and m\/m 2 = (£l/£r) ±;L . For the same 
set of parameters large contributions to the 7 Z channel are possible. For example one can 
reach — 2 for 0l — vr/ 6 , Or ~ vr/10, mi/m 2 — 1.8, m 2 cs 800 GeV and Q ~ 9. This 
region is compatible with S, T and all other constraints. 

A similar analysis can be performed when R w = 2 or larger in Eq. (4.30). In this case 
there are N w pairs of mixing states, with Q = —(N w — l)/2 + Y,... , (N w — l)/2 + Y. For 
each such sector, the mass matrix is 


= ( , (4.40) 

\KQm21 m 2 j 

where is the Clebsch-Gordan coefficient coming from the SU(2) w contraction, deter¬ 
mined by Eq. (B.4), therefore each sector is controlled by the same physical parameters. In 
other words, the two mass eigenvalues, the two mixing angles and the CP-violating phase 
of a given sector determine uni vocally the other sectors too. The corrections to /r 77 and 
p-yZ are obtained summing over the contributions of N w sectors, each being qualitatively 
analog to the case R w = 1 analysed above. Note that, however, one cannot take Q —> 0 to 
recover the SM limit, because there are at least two sectors with different values of Q. Of 
course, the SM is still recovered for small values of the mixing angles. The two fine-tuned 
regions with /U 77 ~ 1 and large pyz are still possible. On the one hand, the different 
sectors can add up to realise AZlw — —2 A/g M - On other hand, for R w = 2 we found a 
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choice of mixing parameters such that ^ vanishes in both the sectors with Q = Y ± 1 /2. 

• Two VLLs coupled to each other, mixing with the SM fermions. Finally, let us 
discuss the possible interactions between ib\ 2 in Eq. (4.30) and the SM leptons. A non-zero 
mixing occurs if and only if 1^1 and/or 1^2 are identified with the states E, L, A or A listed 
in Table 1. There are six such cases, 

R w = 1 : E + L (Y = -1), E + A(Y = 1), 

R w = 2 : L + A (Y = 1/2), A + A(Y = -3/2), A + A G (Y = 3/2), (4.41) 

R w = 3 : A + n {Y = -1), A + fi G (Y = l). 

There are three cases with an additional weak multiplet: a triplet A G = (E, F, G ) ~ 

(1,3,—2), and two quartets 12 = (E c , N, E, F) ~ (1,4,—1/2) and 12 G = (-V, E, F, G) ~ 

(1,4,—3/2), with Q(G) = —3. They do not couple directly to the SM leptons. In these 
three cases the Q = 2 sector couples to the Higgs and therefore may contribute significantly 
to h'y'y and h^/Z. We have shown in section 4.2 that the mixing angles between the SM 
leptons and E, L, A or A must be very small, due to the strong constraints coming from 
the Ztt couplings. With two VLLs the mass matrices become larger, but we expect the 
phenomenology to be qualitatively the same, up to possible fine-tuned cancellations in 
some observable. A crucial effect of the mixing is to make the new leptons decay into SM 
leptons. The components with Q = 2,3 decay more slowly, since their decay chains require 
a virtual exchange of other components of the multiplet. We already reviewed in section 
4.2 the direct bounds on heavy leptons with charge Q = 1,2, decaying promptly into SM 
leptons. We are not aware of any dedicated search for a Q = 3 heavy lepton. 

A detailed analysis of the parameter space is worth only in the context of a specific, 
well-motivated model, and it goes beyond the scope of this paper. The case E + L is 
analyzed in Ref. [74]. The phenomenology of a fourth vector-like family of leptons, L + E + 
N, is studied in detail in Ref. [75]. Here we comment only on the interesting possibility 
to generate the r mass entirely from the mixing with the VLLs, in the limit where the 
SM Yukawa coupling I t lHtr vanishes. There are various ways to induce such coupling 
through mixing, that are illustrated in Fig. 2: 

(o) In the case of E only, one can proceed through a Yukawa coupling connecting I t l to 
Er, followed by two singlet vector-like mass terms. 

(6) Analogously, with L only, one employs two doublet vector-like mass terms and a 
Yukawa connecting Ll to tr. 

(c) In the case of E + L, one can employ a vector-like mass term both for the singlets and 
the doublets, with a Yukawa coupling involving only the new fermions. This case is 
particularly interesting, since it corresponds to the scenario of partial compositeness 
[16] in the r sector: the SM leptons are elementary fields that mix linearly with 
composite VLLs, which couple in turn to a composite Higgs doublet. The SM leptons 
feel EWSB only through the mixing with heavy composite leptons. The Q = 1 mass 
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Figure 2. Different ways to generate the r mass through mixing between the SM leptons and 
VLLs. The dashed lines stand for Higgs vev insertions, the small dots represent a mass mixing 
between a SM lepton and a VLL, and the big dots correspond to the mass of a VLL. The case (c) 
corresponds to the scenario of r partial compositeness. 


matrix and its smallest eigenvalue take the form 
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( 0 tul 0 
0 M l 

\rn E M e 
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(4.42) 


The mixing with the heavy leptons also controls the deviations in the Z couplings, 

2 / \ 2 


$9 


R 


1 ( mE V 

2 \M e ) 


^ 9 TT ~ 


1 / rn L \ 

2 \M l ) 


(4.43) 


As these corrections are bounded by R(Z —> rr), as shown in Eq. (A.27), we find 
that the physical value of m T can be generated for A EE > 2.5, pointing indeed to 
a strong-coupling regime. The phenomenology of r partial compositeness is studied 
e.g. in Ref. [76]. 


(d) Finally, in the case (V'l, V’ 2 ) = (A, E). (L, A) or (A, A), the t mass can be induced by 
three Yukawa couplings, connecting I tE to ^ 2 , tr to ipi, and tpi to ^ 2 , respectively. 
Focusing on (A, A) for definiteness, one finds 


m T 




\av X\v \aav 
\/3Ma \/ 2 M\ \/6 



(4.44) 


For Aaa > 4.5 the physical value of m T can be generated. 


4.4 Vector-like plus Majorana leptons (including higgsinos plus gauginos) 

In this section we consider the interplay between Majorana leptons and VLLs, related 
by one or more Yukawa couplings. If there were no such couplings, the phenomenology 
would reduce to a trivial addition of the effects of Majorana leptons, see section 4.1, and of 
VLLs, see sections 4.2 and 4.3. For reference, the smallest sets of this kind are formed by 
three (four) chiral multiplets: (two copies of) N E or T, E , plus a vector-like pair (El,E e ), 
(Al, Ar) or (A l , A/j). 


• One VLL plus one Majorana lepton. The most general set formed by one Majorana 
lepton interacting with one VLL pair is 

XR ~ (1,-Rtu,0) , ipL,^R ~ 0-1 R w ± 1 ,-1/2) , (4.45) 
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with N v . / 2 + An to avoid the global SU(2) W anomaly. The corresponding Lagrangian is 


- C x ^ = \M x x c r XR + M^uJjr + XiPlHxr + XXrHiI’R + h.c. . ( 4 . 46 ) 

The Majorana mass M x is absent in the case of even SU(2) W representations, N w = 4n: in 
this case there is a conserved ‘new lepton’ number, and all multiplet components combine 
into Dirac fermions. There is a unique physical complex phase, that can be associated to 
M x , choosing M^, A and A real. Each component of the new leptons has a (demi-)integer 
charge Q , therefore it can decay to SM leptons, either through renormalizable interactions 
or higher dimensional operators. For R w = 1 or 3, couplings to the SM leptons can be 
added to the Lagrangian, as we will discuss later. 

In the case ip ~ ( R w ± 1), after EWSB one identifies: one pair of states with charges 
Q = 4z(N w ± l)/2, that combine in a VLL of mass three pairs of states with Q = 
±(N W — l)/2, ±(N W — 3)/2, ..., down to Q = ±1 (±1/2) for N w odd (even); when N w is 
odd, three additional states with <5 = 0. The 3x3 mass matrix in each sector takes the 
form 

( M x kqXv K-qXv\ 


Mq = 


K-qXv 
\ kqXv 


0 

±4 


±4 

o 


(4.47) 


where we chose Q > 0 and a basis with the components of charge Q (-Q) on the left(right)- 
hand side of M.q. The SU{2) W Clebsch-Gordan coefficients k±q are defined by Eq. (B.4), 
and they also include for convenience the factor \j\[2 from the Higgs doublet vev, 


1 lN w + 1 ± 2 Q 

K±Q ~ V2V N ™( N " + 1 ) ' 


(4.48) 


In the case ip ~ (R w — 1), there are three pairs of states with Q = ±(N W — 3)/2, ±(N W — 
5)/2,..., and three states with Q = 0 when N w is odd. The 3x3 mass matrix in each of 
these sectors has the form of Eq. (4.47), but with Clebsch-Gordan coefficients given by 


1 lN w - 1±2 Q 
h±Q ~ v/2 y N W (N W - 1) • 

There are also two pairs of states with Q = ±(1^ — l)/2 and mass matrix 


(4.49) 


( M x_ 

y K-qXv M,p J 


(4.50) 


This is the same structure of Eq. (4.32), that was extensively studied in section 4.3. 

Note that each 3 x 3 or 2 x 2 sector depends on the same five parameters: two masses, 
two Yukawa couplings, and one physical phase. They determine all the mass eigenvalues 
and the mixing matrices. We do not attempt a scan of the parameter space here. The 
corrections to 5, T, /r 77 and /J n z from a sector with the mass matrix of Eq. (4.50) were 
analysed in section 4.3. We expect corrections of the same order from the other sectors. 
Coming to collider searches, for N w ^ 1, 3 there is no mixing with the SM and the lightest 
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mass eigenstate is stable, at least at the renormalizable level. When it is charged, one can 
apply the limits on stable leptons reported in Table 2. For N w odd, the lightest state may 
be neutral, with the typical collider phenomenology of a dark matter candidate. 

The phenomenology is radically different when N w is even. In this case the Majorana 
mass is absent, M x = 0, and the x components are massless before EWSB, therefore one 
mass eigenvalue for each sector is of the order AA v 2 /M^. This implies that all masses are 
bound to the EW scale, whence the situation resembles the one of purely chiral sets of 
new fermions. In particular, one finds a correction to h —> 77 that depends only on Q: 
taking M x = 0 in the mass matrix (4.47) or (4.50), and using the LET approximation of 
Eq. (B.24), we find a CP-even amplitude q — 8Q 2 /3 for a 3 x 3 sector, and the same 
for a 2 x 2 sector. For comparison, when N w is odd the Majorana mass is allowed and one 
finds 


< 77 , 3 x 3 




8 Xv Xv / 9 

- Q (k z c 

3M X M^ 1 ~ Q 


+ k Q ) > 


< 77 , 2 x 2 




8 Xv Xv 
3 M x 


QV 


Q > 


(4.51) 


where we made the approximation Xv, Xv -C M x , The correction grows as the x ~ ip 
mixing parameters of the type A v/M increase, on the other hand S and T generally put a 
significant upper bound on these parameters. Obviously, in each model the total Anew is 
the sum over all the sectors with different Q. 

Let us say a few more words on the cases where x and/or if} mix with the SM leptons. 


• N + L : for R w = 1, the Majorana fermion is a sterile neutrino N, the vector-like 
fermion is a lepton doublet L, and the Lagrangian in Eq. (4.46) is extended to include 
XnIlHNr + XlLlHtr + Ii.c.. The full parameter space includes two real masses, four 
reals Yukawa couplings and two physical phases, that can be associated e.g. to A and 
A. The 2x2 mass matrix in the Q = ±1 sector is given by Eq. (4.13): the mixing with 
the SM is small due to the strong constraint from the Ztt couplings, implying small 
deviations in and ii-,z■ The 4x4 mass matrix in the Q = 0 sector is obtained 
by adding to Eq. (4.47) a first row and a first column of the form (0, X^v/V2 ,0, 0). 
There are some simple limiting cases. If A, A —> 0, the matrix reduces to two diagonal 
blocks, as N and L decouple from each other, and the phenomenology reduces to the 
one of the previous sections. If AIl,n/v are much larger than the Yukawa couplings, 
all the mixing angles are small and the smallness of the light neutrino mass follows 
from the usual seesaw mechanism, m u — X 2 N v 2 /{2 Mjv). Still, m u can be small even 
in the presence of large Yukawa couplings. In particular, for Xv,Xv M NtL , large 
mixing angles are possible, and, correspondingly, such intricate neutral sector may 
induce significant corrections to the S and T parameters. 


• For R w = 3, the Majorana fermion is the triplet E and the vector-like fermion is the 
doublet L or the quartet fl. 

E + L : the Lagrangian in Eq. (4.46) is extended to include X-^IlHTir + A lLlHtr + 
h.c.. Up to different Clebsch-Gordan coefficients, the 4x4 neutral sector is the same 
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as in the N + L case. The mass matrix in the charged sector is 


Mi = 



0 ^ 
Ms yj\.Xv 
\f\xv M l j 


(4.52) 


The Ztt couplings constrain both mixing parameters Asn/Ms and A lv/Ml to be 
small, as explained in sections 4.1 and 4.2, respectively. When the mixing with the 
SM is neglected, one is left with a special case of Eq. (4.50), which corresponds 
to the chargino mass matrix in supersymmetry. Note that, in the limit where A T 
vanishes, there is still a contribution to the r mass, m T — (AsA^Au 3 )/(3\/2MsM£), as 
illustrated in Fig. 2(d). Despite the constraint on the mixing from the Ztt couplings, 
one can accommodate the correct size of m T for A > 3, see the discussion below 
Eqs. (4.42) and (4.44). 

E + D : the Lagrangian in Eq. (4.46) is extended to include Xy,IlHY> r + h.c., as the 
quartet does not mix with the SM. The 4x4 neutral sector has the same structure 
as in the N + L and E + L cases. The Q = 1 sector also has a 4 x 4 mass matrix, that 
is obtained from Ecp (4.47) by adding a first row (A T u/\/2, Ae/\/ 3, 0,0), and a first 
column (A T u/\/2,0, 0, 0). In addition, there is a Q = 2 state with mass Mq. There 
are no significant phenomenological novelties, as the effects of the SM mixing with E 
and of the E mixing with D mixing do not interfere significantly. 


• One VLL plus two Majorana leptons. Let us come to sets of two Majorana leptons 
both interacting with one vector-like pair. One obvious possibility is to take two copies of 
the same Majorana lepton, that is, to replace xr i n Eq. (4.45) with XiRi * = 1,2, with the 
obvious doubling of each coupling involving \ i n the Lagrangian. Note that N w can be 
arbitrary and, for even N w , the Majorana mass terms are forbidden but a Dirac mass term 
M xX I rX2 r is allowed. In all other respects, the mass matrix structures and the inherent 
phenomenology are a straightforward generalisation of those discussed above. 

The second and last possibility to couple two Majorana leptons to one VLL is provided 
by the set 


XlR ~ (1, Rw , o) , V’L, IpR ~ (1, Rw + 1, -1/2) , 


X2R ~ (1, Rw T 2,0) , (4.53) 


with N w necessarily odd, and Lagrangian 

- = M^ L ip R + Y?i=i \ M XiX c iR XiR + Xi'i/jLHXRi + XiX c Ri H^ R + h.c. . (4.54) 


There are two pairs of states with Q = ±(N w J r l)/2, with a mass matrix given by Eq. (4.50) 
with M x , A, A —> M X2 , A 2 , A 2 . In addition, there are four pairs of states with Q = ±(N W — 
l)/2, ±(1V W — 3)/2,..., ±1, and four states with Q = 0. The 4x4 mass matrix in each 
such sector takes the form 


Mq 


( M Xl 0 ki,q\iv ki -qAiu^ 

0 M X2 R2,qX 2 v K2-q\ 2 V 

ki-qXiv k 2 -qX 2 v 0 
\ ki,qXiv k 2 ,qX 2 v M,p 0 / 


(4.55) 
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with k\,±q given by Eq. (4.48), and K 2 ,±q given by Eq. (4.49) with N w —»• N w + 2. As in 
Eq. (4.51), one can estimate the contribution of this mass matrix to the /177 coupling, by 
taking the LET approximation, 

477,4x4 ^ 

"W1X2 ^,Q 

Let us discuss the minimal cases R w = 1 and R w = 3, that admit a mixing with the 
SM leptons. 

• N + L + E : for R w = 1, the new fermions have the gauge quantum numbers of the 
bino, the higgsinos and the wino in supersymmetry. Thus, the mass matrices (4.50) 
and (4.55) are a generalisation of the chargino ( Q = ±1) and neutralino (Q = 0) 
mass matrices, respectively (for a review see Ref. [77]). Supersymmetry restricts 
the Yukawa couplings to Ai^/Ai^ = —tan (3, A 1 /A 2 = A 1 /A 2 = —tan 9 W /V3, and 
Ai = — cos f3g' . The effect of charginos and neutralinos on the Higgs boson couplings 
is analysed e.g. in Refs. [78-81]. In particular, the chargino loop contributing to /177 
and h'yZ is controlled by the weak coupling g, and it is typically subleading compared 
to the SM top quark loop. Without the supersymmetry constraints, the most general 
chargino mass matrix has the structure of Eq. (4.32), therefore one can apply the 
results of section 4.3 for the Higgs decay amplitudes into 77 and 7 Z. 

In the absence of (A-parity conserving) supersymmetry, not only the four Yukawa 
couplings A 12 and Aq 2 are unconstrained, but in addition a mixing with the SM lep¬ 
tons is allowed: one should add to the Lagrangian in Eq. (4.54) the terms XnIlHNr + 
\'zlLHYjR+\] J LLHTii+h.c.. The Q = ±1 mass matrix becomes the one in Eq. (4.52), 
and the Q = 0 mass matrix becomes 5x5, and it is obtained from Eq. (4.55) by 
adding a first row and column of the form (0, Aatw/\/2, \t,v/VQ, 0, 0). Therefore, one 
can observe the phenomenological effects of N, E and L individually, as analysed in 
sections 4.1 and 4.2, as well as their interplay, already described above for (N + L ) 
and (E + L). As usual, the mixing with the SM leptons is typically constrained to 
be small by the smallness of m u and by the Ztt couplings, thus the modifications 
to the hw and hrr couplings are suppressed. However, even a very small mixing 
with the SM offers decay modes to the heavy fermions, such that none is stable. A 
dedicated analysis of the full parameter space would be interesting, to characterise 
quantitatively the correlations among the different observables, and especially the 
deviations from the supersymmetric limit. 


A iV A iV 

M Xi M^j, 


Q 2 «-Q + 4q) 


(4.56) 


• E + fI + 5: for R w = 3, the new fermions are a Majorana triplet, a vector-like quartet, 
and a Majorana quintuplet S. There is a 2 x 2 sector with Q = ±2 given by Eq. (4.50). 
As the triplet mixes with the SM through A y,IlHY>r + h.c., there is a 5 x 5 sector 
with Q = ±1, that is obtained by adding to the matrix in Eq. (4.55) a first row 
(A t u/\/2, X^v/y/3, 0, 0,0) and a first column (A T u/\/2, 0, 0, 0,0). This large number 
of charged states with potentially large mixing can give a significant correction to 
// 77 and n<yz- For concreteness, neglecting the mixing with the SM and using the 
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LET approximation, we find 


All ^ 

3 


1 Ain Aye X2V X2V 
3 Ms Mq Mg Mq 


(4.57) 


to be compared with the SM top contribution, A/ 7 — 16/9. One should take into 
account the constraints (in particular S and T) on the mixing parameters ~ A v/M. 
The neutral sector has also a 5 x 5 mass matrix, obtained by adding a first column 
and row (0, X^v/y/6, 0, 0,0) to to the matrix in Eq. (4.55). As usual, the vanishing 
neutrino mass requires As^/Ms to be very small. 


5 Phenomenology of non-chiral quarks 


In this section we discuss new coloured fermions that either form vector-like pairs, or 
admit a Majorana mass term. We will dub them ‘quarks’ even when they are not in the 
fundamental representation of SU( 3) c . 


5.1 One vector-like quark 

There are seven possible VLQs that mix with the SM quarks, as listed in Table 1. They 
have been extensively studied under various respects in the literature (see e.g. Refs. [22, 82- 
87]). Here we describe in a compact, systematic way the leading order constraints coming 
from EW precisions tests, direct searches at colliders, and Higgs couplings. As usual we 
restrict ourselves to mixing with the third family. In the top (bottom) sector, a mixing 
appears whenever the VLQ contains a component T ( B ) with the same charge as t ( b ). 
The components of each multiplet are displayed in Eq. (2.6). In the case of weak singlets 
or triplets, the SM Lagrangian is extended by 

- = X^qZH{H)ifi R + M^ L 'i/; R + h.c. , for = T,Xq(B,Yq) , (5.1) 


and, in the case of weak doublets, by 

-£* = XfifoH(H)t R + X b ^T L H{H)b R + + h.c. , for i/j = Q(X T , Y B ) , (5.2) 

with the further restriction A y t = Ay = 0. The structure of the top (bottom) sector mass 
matrix is very close to the charged lepton one in the case of one VLL, therefore we will 
frequently refer to section 4.2. In the top sector one has 


•* = T ' x *r»’ 

with Clebsch-Gordan coefficients k-tq x T = 
rotation to the mass basis is parametrized as 
Eq. (4.14). In the bottom sector one has 


M,= 


^73 u 

,44^ 


, ^ 


1 , i\ Xq = yM/3 and k 
M t = U L diag(m tl m t ')U 


t 

yQ 

t 

w 


= Q,x t , 

= 4273 . 

in analogy 


(5.3) 

The 

with 


M b 


(X b f 2 4A^\ 

V0 ) 


^ = B,Xq,Yq , Mb 


(M 

V44t2 


^ = Q,Y b , (5.4) 
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with n b B qy b = 1, K b x Q = -\/2/3 and Ky Q = \/l/3, and one can write Mb = U B diag{mb, rriy)U' R 
In all cases except Q , the vector-like mass M,j, and the three independent Yukawa 
couplings can be taken to be real. In the case of Q, there are four Yukawa couplings and 
one complex phase 4> is physical. In full generality, one can choose A t) b and X*q real, and 
add a matrix P ,A = diag(e u 1) on the left of Mb- Then, Eq. (A.4) and Eq. (B.5) show 
that (j) appears in the W couplings to the fermion mass eigenstates, while the Z and h 
couplings are independent from </>. It is very difficult to observe such CP-violating effect 
in the charged current, since it vanishes for —>• 0, that is, it is always suppressed by the 

small ratio mb/v. Coming to the mixing angles, the left- and right-hand ones are related 
as in the case of VLLs, see Eq. (4.15): in the case of weak doublet VLQs one finds 

tan 0 B = —-tan0R < tan0R , tan 0 B = — -t&nO B <C tan Or , (5-5) 

my my 


while for singlet or triplet VLQs, the same relations hold with L <->• R. In the following 
we will drop the subscripts L , R and denote 9^ (6-^) the largest mixing angle in the top 
(bottom) sector, for any given VLQ if). For the multiplets with both the T and B com¬ 
ponents, the mass eigenvalues and the mixing angles in the top and bottom sectors are 
strictly related, 

{ 2 2 2/2 2 \ ~ 2/2 2 \ 
m t , - m b , = - m t ) - s^{mp - m b ) 

- mj)r^ = s^(m^, - m 2 b ) 

where rx Q = \/2, t Yq = l/y/2, and rg = Xq/Xq. Therefore, one can determine the 
bottom sector parameters, my and 9^,, as a function of the top sector one, my and 0^, or 
vice versa. In the case = Q, there is the additional freedom of the choice of rg. Note 
that the custodial symmetry is preserved in the Q sector if rg = 1, see Eq. (5.2). The mass 
splitting among the heavy quarks is controlled (at tree level) by the mixing with the SM. 
The mass ordering is determined as 


T : Mt < my , B : Mg < m-y , 

X T : m x = M Xt < my , Y B : my = M Yb < my , (5.7) 

Xq : m x = M Xq < my < my , Yq : m Y = M Yq < my < my , 


where we took implicitly into account the experimental upper bounds on the mixing and 
on mt/my, when needed to establish the ordering. In the case of Q, for rg = 1 one finds 
Mg < my < my, but the ordering between b' and t' can change for different values of rg. 

Masses and mixing angles are constrained by the perturbativity of the Yukawa cou¬ 
plings, 


\Kl >| - 


V2 


my 

t s i > 

v 


<C 47t 


\ b 


y/2 _ m-y 



b~ s tl> 

4 v 


■C 47t . 


(5.8) 


Note that we do not impose a stronger upper bound such as Air/y/W c , for reasons discussed 
in Appendix B.l. The perturbativity of the SM couplings A* and X b is guaranteed a 
fortiori. For definiteness, in Figures 3 - 6 we delimit with a black dotted line the region 
of parameters where at least one Yukawa coupling becomes larger than 27r. In the case of 
Q both inequalities must be satisfied at the same time, therefore a large departure from 
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Figure 3. Constraints on the weak singlet VLQs B (left panel) and T (right panel), as a function 
of the mass of b' (F) and of the mixing angle between by, and b' L (tL and t' L ). The region above 
the dotted black line is excluded by perturbativity. The blue-shaded region is excluded by the Zbb 
couplings. The (light) yellow-shaded region is excluded by the S and T parameters at (68%) 99% 
C.L.. The green-shaded region is just the intersection of the previous two. The grey-shaded region 
is excluded by the collider searches summarised in Table 3. The region above the solid black line 
is excluded by a rough global fit of the Higgs couplings at 99% C.L.. The dashed blue (dotted red) 
lines correspond to a few relevant values of the signal strength fijz (g-yy)- 


rQ = 1 leads to a stronger constraint, as illustrated by the comparison of the left and right 
panels of Fig. 6. 

Important constraints come from the Z couplings to quarks, that are affected by the 
mixing as shown in Ecp (A.6). The tree-level deviations with respect to the SM are given 
by 


^ g Zbb 


= 4 


2+ T b 


^9zhh ~ Sr t b 


* Zbb 


^9ztt ~ S L 


- 2 +n 


5 9ztt ~ s R T t 1 
(5.9) 


where Tg (Tj,) is the weak isospin of the B (T) component of the VLQ under investigation. 
The Zbb couplings are measured less precisely than their leptonic analog, the Zt + t~ cou¬ 
plings, but nonetheless they are strongly constrained, especially for 6^. The top couplings 
to the Z boson are poorly constrained directly, however they also contribute at one loop to 
Zbb. These constraints are summarised in appendix A.3 and they exclude the blue-shaded 
region in Figures 3-6. Deviations in Zb^bL are present at tree level in the case of B, Yb , Yq 
and Xq. However, in the case of the doublet Yb, the deviation is suppressed by ( mb/m' b ) 2 
and the most important correction is the one to Zb^bR. In the case of T and Xt, there 
is no bottom partner and the deviation to Zb^L is induced at one loop mostly through 
bfJz t t , leading to a relatively weak constraint. Finally, in the case of Q, Sg ^ is generated 
at tree level and 5g^ b at one loop, the strongest constraint coming from the right-handed 
coupling. 
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m b - (GeV) m,< (GeV) 

Figure 4. Constraints on the weak doublet VLQs Yb (Y = —5/6, left panel) and Xt (Y = +7/6, 
right panel) as a function of the mass of b' (t') and of the mixing angle between bn and b' R {tn and 
t' R ). The notation is the same as in Fig. 3. 


The VLQ couplings to the EW gauge bosons are also constrained by the EW precision 
parameters S and T, whose expressions are provided in appendix A.2. Note that, in 
contrast with the tau and bottom sectors, in the top sector contributions proportional to 
powers of mt/mz are not suppressed. In Figs. 3 - 6 we display in (light) yellow the region 
corresponding to the ( 68 %) 99% C.L. ellipse of Fig. 10. Since S and T are proportional to 
the mixing between the SM quarks and the VLQ, one typically observes an upper bound 
s. 0 , < 0.05 — 0.20, depending on the VLQ under consideration. Note that this bound 

relaxes as the heavy quark mass decreases, because S and T eventually vanish in the limit 
mt> —>■ mt or my -+ m^. Note also that a cancellation is possible among relatively large 
contributions to S and T, such that large mixing angles may be allowed in a fine-tuned 
region of parameters. This is especially relevant in the case of Xt, because such region is 
not excluded by other constraints. Indeed, we find 
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(5.10) 

where we dropped terms subleading in sx T and mt/my. As the logarithm is large and 
negative, a cancellation is possible in the T parameter even for large mixing: this explains 
the allowed strip in Fig. 4, that reaches sx T — 0.5. A comment is in order for the case 
of Q : one would expect a milder constraint from T when the two Yukawa couplings Aq 6 
respect the custodial symmetry, i.e. when tq = 1. However, even in this case there is 
an important deviation from the SM, because the residual custodial-breaking parameter, 
(A* — A b), differs from the SM one, y/2(mt — mb)/v, as soon as the mixing is non-zero. 

Let us now turn to the direct searches of VLQs at colliders. As they are coloured, it is 
easier to produce them at the LHC, relatively to VLLs. Below ~ 1 TeV they are dominantly 
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Figure 5. Constraints on the weak triplet VLQs Yq {Y = —1/3, left panel) and Xq (Y = +2/3, 
right panel) as a function of the mass of t' and of the mixing angle between and t' L . The notation 
is the same as in Fig. 3. 




m r (GeV) m r (GeV) 

Figure 6. Constraints on the weak doublet VLQ Q as a function of the mass of t' and of the 
mixing angle between t R and t' R . The left panel corresponds to Xq = Xq {tq = 1), and the right 
one to 4 Aq = Xq (tq = 4). The notation is the same as in Fig. 3. 


produced in pairs through strong interactions, while for higher masses single production 
by EW interactions can become dominant [83, 88]. The pair production mechanism, that 
dominates in the mass range probed at the 8 TeV LHC, is independent from the VLQ (all 
are colour triplets) and from the mixing parameters. The ATLAS and CMS searches focus 
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on the following decay channels for the heavy quark mass eigenstates: 


X->tW + , t' -*tZ,th,bW + , b' —»• bZ, bh, tW~ , Y -+ bW~ . (5.11) 


Note that t' and b' can decay via neutral current at leading order, owing to their vector¬ 
like nature. Decays into another heavy quark coming from the same multiplet, such as 
t' —y XW~, are kinematically suppressed; decays through loops may also be possible, but 
they are typically negligible [84], Here we will disregard these sub-leading channels, and 
assume that the decay channels in Eq. (5.11) have unit branching ratio. In addition, the 
decays are assumed to be prompt, that is the case whenever the mixing angles are large 
enough to have an observable effect on the Higgs couplings. 

The relative branching ratios of t! and b' in the three decay channels depend mostly 
on the weak isospin of the VLQ and on the mixing angles. Indeed, since the heavy quarks 
are already constrained to be heavier than a few hundred GeVs, in good approximation 
one can neglect the final state masses and find 


r(f' -+ ht) ~ (\y tt '\ 2 + \y tt '\ 2 ) , 

r(f' ->• zt) ~ 327r ^| (K c fM 2 +1( c iM 2 ) > 

r(t ' wb) * 3^% (l^l 2 + 1( c SVl 2 ) 


(5.12) 


and the same for t! o b' and t o b. Here the Higgs couplings are defined by Eq. (B.5), 
and the Z and W couplings by Eq. (A.4). By a straightforward computation, one finds 
Br (t' —> ht) — Br {t’ —> Zt) and 


T : Br (t' —> Zt) — 
X T : Br (tJ -> Zt) ~ 
Q : Br (t' -> Zt) ~ 

Aq : Br {f ->• Zt) ~ 
Eq : Br (f Zt) ~ 
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(5.13) 

As before, we neglected the SM masses and, therefore, the subdominant mixing angles in 
Eq. (5.5). Note however that some branching ratios are proportional to the SM masses at 
leading order, for example in the case of Xt one finds Hr (t 1 —> Wb) ~ m 2 /(c\ T m 2 ,). In 
the cases where both t' and b' are present, we used the relation tan 0^ ~ r^t&nQ^, that 
follows from Eq. (5.6) if one neglects mb and rrq. The b' branching ratios are obtained from 
Eq. (5.13) by the replacements T —> H, Xt —5- Yr, Xq -h- Yq , t' — > b ', t -f-)- b, rQ —> 1/rg 
and Sip —>• 4 for each ?/>. 

The experimental lower bounds on the mass of t! and b' are presented as a function 
of two independent branching ratios [89, 90]. We choose, in the plane of branching ratios, 
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heavy quark 

branching ratios 

multiplets 

mass bound 

V (Q = 5/3) 

Br wt = 1 

X T , X q 

rrix > 840 GeV [91] 

t' (Q = 2/3) 

Br W b = \, Br zt = Br ht = \ 

T, X Q , Q(r Q = 1) 

rn t > > 800 GeV [89] 

Br Wb = 0, Br Z t = Br ht = \ 

Xt, Yq, Q( r Q < 1) 

rnf > 855 GeV [89] 

Br Wb = 1, Br zt = Br ht = 0 

Q{tq » 1), Q + Y b 

rn t i > 920 GeV [96] 

Br Wb = Brzt = 0, Br ht = 1 

X t + Q 

rn t > > 950 GeV [89] 

Br Wb = Br ht = 0, Br Z t = 1 

X t + Q 

m t , > 800 GeV [89] 

Br Wb + Br Z t + Br ht = 1 

T, X T , Q, Y q , X q 

rn tJ > 720 GeV [96] 

b' (Q = -1/3) 

Br wt = Br Zb = Br hb = \ 

B, Yq, Q(t q = 1) 

m b , > 735 GeV [89] 

Br wt = 0, Br Zb = Br hb = \ 

Yb, Xq Q(r 2 Q » 1) 

my > 755 GeV [92] 

Br wt = 1, Br Zb = Br hb = 0 

Q(vq <C 1), Xt + Q 

my > 810 GeV [91] 

Br wt = Br Zb = 0, Br hb = 1 

Q + Y b 

m b ' > 846 GeV [94] 

Br wt = Br hb = 0, Br Zb = 1 

Q + Y b 

my > 775 GeV [89] 

Br wt + Br Zb + Br hb = 1 

B, Y b , Q, Yq, Xq 

my > 582 GeV [95] 

Y(Q = -4/3) 

Br Wb = 1 

Yb, Yq 

m Y > 770 GeV [89] 


Table 3. Lower bounds at 95 % C.L. on the heavy quark masses mx,mt',rn b ’ and my. The 
experimental searches assume pair production via strong interactions and prompt decays in the 
indicated channels. In the second column we specify the assumption on the heavy quark-decay 
branching ratios. Here Brzt stands for Br (t' —> Zt ), and so forth. In the third column we list 
the VLQ multiplets that correspond to those branching ratios, in the small mixing approximation. 
Here “Xt + Q” and “ Q + Yb ” refer to pairs of VLQs with a custodial symmetry, that are discussed 
in section 5.2. 

their approximate values for the VLQ under consideration. To this purpose, we take the 
limit Sip —y 0 in Eq. (5.13) (3/. —> 0 in the case of b'), because the collider bound is relevant 
at small mixing angles, see Figs. 3-6). The only exception is Xt , where large mixing is 
possible, but in this case the strongest collider bound is the one on the component X. The 
lower bounds on each heavy quark mass are collected in Table 3, and vary between ~ 600 
and 900 GeV [89, 91-95]. The region excluded at 95 % C.L. is shaded in grey in Figs. 3-6. 
A detailed analysis of the lower bound on m' t is presented in Ref. [23] for the case of the 
VLQ T, taking also into account indirect constraints from R-physics observables. 

Let us now discuss the corrections induced by the VLQ on the Higgs boson couplings. 
The couplings of t,t',b and b' to the Higgs have the same form as those of t and t' in 
Eq. (4.19), with the obvious replacement of masses and mixing angles. The heavy quarks 
X and Y do not couple to the Higgs. The Higgs signal strengths at the LHC pt a , defined in 
Eq. (B.40), are the product of three factors: Higgs production rate, partial decay rate and 
lifetime. While new leptons affect significantly the partial decay rate only, new quarks can 
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modify substantially each factor. In particular, the Higgs production via gluon fusion is 
sensitive to a VLQ, and the mixing in the bottom sector can change significantly the total 
Higgs width IV Let us remind that, as discussed below Eq. (5.4), the Higgs couplings are 
CP-conserving for any VLQ. 

When the VLQ contains a B component, the Higgs width into bb is modified, with 
respect to the SM, by a factor R^ = (1 — s^) 2 that, in the light of previously discussed 
constraints, can be as small as ~ 0.9. Since in the SM h —> bb is the dominant decay channel, 
this correction enhances g a for all other decay channels, through the factor Tf t M /Th- The 
Higgs production via gluon fusion is modified by a factor 


R 


J 99 


A 3 SM + \ s % [^1/2 ( r i') - M/2 ( r t)] + 1^; [^ 1/2 ( T b') - Ai/ 2 (Tb)\ 

I AliA 2 


(5.14) 


where the form factors are defined in appendix B.2.1. The effect of the top sector is 
qualitatively different from the bottom one: given the collider lower bound on rrif/ y , their 
form factors are very close to the asymptotic value, H 1 / 2 (0) = 4/3. While the t loop is 
also close to this value, the b has a small mass and a suppressed form factor: H 1 / 2 (^') — 
A 1/2 (r t ) ~ —0.04 and A±/ 2 (tv) — ^4 1 / 2 (u) — +1.41 Therefore, when a b' is present (for i/j = 
B,Yb,Q, Xq,Yq), its effect dominates and the interference with the SM is constructive. 
An exception is possible for -0 = Q, where sq/sq <C 1 for vq <C 1, see Eq. (5.6). In the 
latter case, and when only a t! is present (for i/j = T,Xt), there is a slight destructive 
interference with the SM. The tth production mode is also modified respect to the SM in 
the presence of t — t! mixing, with a cross-section reduced by a factor . 

In the diphoton channel 




AJm + [^ 1 / 2 ( T t') ~ A 1 / 2 (T t )\ + \Ai/ 2 (tv) - Ai/ 2 (n)\ 


14 


77 | 2 
SM I 


(5.15) 


Here the SM amplitude is negative, therefore the interference pattern is reversed with 
respect to R gg . A few different values of the signal strength // 77 are shown in Figs. 3-6 
by dotted red lines. Once the other constraints are taken into account, one finds at most 
<V 77 ~ 0.3. Finally, for the Higgs decay into a photon and a Z, the new physics amplitude 
writes 


a=t,b 


A? m+ , - AA s Z m ~ ^ —^A-^5g^ a Ai/ 2 {r a , A q ) + s^ a [ga'a'Ai/ 2 (r a >, \ a >) - g^ a A 1 / 2 {r a , A a )] 

(5.16) 


m a + m a ' y , , x \ 

+ ViO+V 0 /- 5oa , '^l/2 vc a' i ^a' i ha i ) r i 

yjT\'l GL Tfl GL i J 


where s ^ = s^, = Sy, and we neglected the terms proportional to the form factor 

B]_/ 2 i that are subdominant. Note also that there is no CP-violating amplitude, as both 
the h and Z couplings respect CP, see discussion below Eq. (5.4). The vector Z couplings 
are obtained from Eq. (A.6): one Ends 5g^ a = (g^ a - gX’a ) ~ s^ a and g v aOL , ~ s^ a , so 
that the new physics amplitude of Eq. (5.16) is of the order of the mixing squared. The 
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interference with the SM may be constructive or destructive depending on the sign of the 
Z couplings. A few different values of the signal strength p,^z are shown in Figs. 3 - 6 by 
dashed blue lines. When a b' is present, both p^z and /r 77 receive a similar correction, 
dominated by the increase of R gg and Ff t AI /T^. A significant correction is possible for 
Yb, with both 5p^z and J/r 77 as large as ~ 0.3. On the other hand, in the case of T or 
Xt the corrections to the two channels are significantly different, because the small factor 
.Ai/ 2 (tv) — Aij 2 (r t) in Eq. (5.15) suppresses the correction to L 77 . In particular, for Xt 
one can have 5p,jz ~ 0.1 with <5/r 77 ~ 0 . 01 . 

In our analysis we computed the relevant signal strengths p a for a = bb, 77 , 7 Z, WW 
and ZZ, taking of course into account the corrections to R gg and F^. We compared these 
predictions with the allowed experimental ranges given in Table 4. By simply computing 
a x 2 f° r these five channels, we determined the region of parameters disfavoured at 99% 
C.L., that is delimited by the solid black line in Figs. 3-6. 

5.2 Two vector-like quarks (including b and t compositeness) 

• Two VLQs not coupled to each other. Let us consider first a pair of VLQs (ip,4> f ) 
that do not couple to each other via a Yukawa coupling. In this case ?/> and ?// must be 
identified with one of the seven VLQs in table 1, and the phenomenological effects are, 
in most respects, a trivial addition of those of each VLQ separately, already discussed in 
section 5.1. 

A noticeable exception occurs when the Yukawa couplings and the vector-like masses of 
-0 and i/j 1 respect an additional SU(2) R global symmetry, that provides custodial protection 
for the EW gauge boson couplings: when the parameters approach this custodial limit, the 
constraints from EW precision tests drastically relax with respect to the case of a unique 
VLQ. There are four pairs that may form a doublet under SU (2)#: the weak singlets (T, B ), 
the weak doublets (Xt, Q ) or ( Q , Yb), and the weak triplets ( Xq, Yq). For illustration, we 
will concentrate on the case of doublets. 

The two VLQs transform as bi-doublets under a custodial SU ( 2 )l x SU(2) R symmetry, 
as long as their Yukawa couplings to the SM fermions and their vector-like masses are equal, 

~£(x t ,q) = ( x t Q)l ( ft ) t R + (Xt Q)l ( I J + h.c. , 

> / f\ R (5 ' 17) 

-£(Q,y b ) = ^TQYbVl b R + M^(QYb)l ^j + h.c. . 

These are the smallest sets of top and bottom quark custodians, respectively [70, 97, 98]. 
Note that, in this custodial limit, the additional coupling Q L Hb R (Q L Ht R ) must vanish 
in the top (bottom) case. Therefore, a mixing occurs only in the top (bottom) sector, and 
there are no deviations in the bottom (top) couplings, despite the presence of a b' (t') in 
the spectrum. The analysis is analogous to the case of r custodians, discussed in section 
4.3. For example, in the top case the linear combination t" = (T (Xr r) _ T^)/\/2 does not 
couple to the Higgs and therefore it does not mix, while the orthogonal combination mixes 
with the SM top quark as in Eq. (5.3), to form the mass eigenstates t' and t. The mass 
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Figure 7. Constraints on the pairs of VLQs ( Q , Yg) (left panel) and (Xt, Q) (right panel), in the 
custodial limit (equal vector-like masses and Yukawa couplings), as a function of the mass of b' (t') 
and of the mixing angle between b g and b' R (t R and t' R ). The notation is the same as in Fig. 3. 


spectrum is 


Mj, 

mt" = mb' = rrix = Mj> < mf — - 

(X T ,Q) 

Mj, 

mb" = mf = my = Mjj < mb' — —— 

CR 

(Q,y b ) 


(5.18) 


Due to the custodial symmetry, the values of the heavy quark branching ratios into SM 
particles differ from the case of a single VLQ, discussed in section 5.1. We assume that 
the decays to another heavy quark are kinematically suppressed, because of the small mass 
splitting in Eq. (5.18), and once again we neglect the SM masses in the final state, as well 
as the t — t' (b — b') left-handed mixing angle, that is suppressed by mt/m# (mb/my). 
Consider for example the top case. The decays t'/t" —> Wb are suppressed as in Eq. (5.13) 
(here tq = 0). In addition t' —> Zt vanishes because the Q = 2/2, components of Xt and Q 
have opposite weak isospin, and t" —> ht vanishes because t" does not couple to the Higgs. 
Similar arguments hold in the bottom case. In summary one finds 


BR(t" -> Zt) ~ Br(b' -»• W~t) ~ Br(X -> W+t) ~ Br(t' ht) ~ 1 (X T , Q) , 

BR(b" -f Zb) ~ Br(t' -> W+b) ~ Br(Y -> W~b) ~ Br(b' -t hb) ~ 1 (Q, Y B ) ■ 

(5.19) 

It is amusing that, in these two models, there is one heavy quark decaying exclusively in 
each of the possible decay channels listed in Eq. (5.11). The experimental lower bounds on 
these heavy quark masses can be read off Table 3. 

The custodial symmetry protects the Zbb couplings: in the top case 5gj~ = 0 and 
the small dgj^ = (mt / M^,) 2 s 2 R / 2 contributes to Zbb only at one loop; in the bottom case 
5gj^ = 0 and 5gj^ = (mb/M^) 2 'sj i /2 is very suppressed by the bottom mass. Thus, in this 
custodial limit large mixing angles are not excluded, as shown in Fig. 7. Indeed, one can 
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see that the constraint from the other EW precision parameters is significantly relaxed 
too, as T receives a small correction only, from the difference (At — Aft), that is not SM-like 
because of the mixing, while S acquires a positive correction that remains in the ellipse 
unless the mixing is very large. Note that in the bottom sector the T parameter is al¬ 
most independent from my, since the smallness of mj, ensures Aft — ■ On the other 

hand, when the mixing occurs in the top sector, the coupling A t and consequently the T 
parameter strongly depend on mf. The dominant constraints at small and large heavy 
quark masses come from the direct collider searches and from perturbativity, respectively. 
In the top case (right panel of Fig. 7), the mixing is not constrained by the fit of the Higgs 
coupling, as the bottom sector is SM-like. As a consequence, for 1.5 TeV < m t ' < 2 TeV 
the mixing can be as large as sr ~ 0.8. The 77 Higgs decay channel can be suppressed 
at most by <5 // 77 ~ —0.03, while the corrections to the 7 Z channel may be larger, up to 
5fi<yz — +0.13. In the bottom case (left panel of Fig. 7), the mixing in the bottom sector 
enhances all the other Higgs channels, as R^ = (1 — S|j ) 2 < 1. This leads to an upper 
bound sr < 0.35. Significant corrections as large as <5 // 77 ~ — 0.6 are possible. Note 

that the Higgs signal strengths in Fig. 7 are similar to those with Yr or Xr only, shown 
in Fig. 4. The difference is that the region allowed by EW precision tests largely inflated 
here, thanks to the custodial symmetry. 

• Two VLQs coupled to each other, not mixing with the SM fermions. Let us 

move to the case of two VLQs coupled to each other via Yukawa interactions. Their chiral 
components transform as 

VhZo i>iR ~ (Rc,Rw, Y) , ^2L,^2i? ~ (R c ,Rw + 1,H + -) , (5.20) 

with R c 1. The corresponding Lagrangian is the same as in the case R c = 1 (two 
VLLs), and it is given in Eq. (4.31). We bar the special case R c = R c , Y = 0 and R w odd 
(Y +1/2 = 0 and R w even), that allows for a Majorana mass terms for 7/7 (V’ 2 ) and will be 
discussed in section 5.3. We also bar mixing with the SM quarks, that will be discussed at 
the end of the section. The effect of two VLQs on /r 77 was discussed in detail in Ref. [99]. 

The number of mass eigenstates with a given electric charge Q and the structure of 
their mass matrices are the same as in the case of two VLLs, see Eqs. (4.32) and (4.40). 
Therefore, there are five physical parameters: two masses 7711 , 2 , two mixing angles 6l,r 
and one phase <p, defined by Eq. (4.33). The analysis of the parameter space proceeds 
exactly as in section 4.3 and will not be repeated here, however the phenomenology is 
strongly modified as the colour representation R c is non-trivial. The main differences are 
the following: 

• The VLQs are pair-produced via strong interactions and, in the absence of mixing 
with the SM, the lightest state is stable and hadronises. The direct collider bounds 
on these particle masses are above one TeV, as we already described in some more 
detail in section 3. 

• The contributions of the VLQs to the S and T parameters, as well as to the Higgs 
decay amplitudes into 77 and 7 Z, have the same form as in Eqs. (4.35) to (4.39), 
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Figure 8. In the left (right) panel, we show the constraints on the pair of VLQs B and Q (T and Q), 
in the limit of vanishing mixing with the SM quarks, degenerate mass eigenvalues my = mb" = m \/3 
(m*/ = mt" = TO 2 / 3 ) and no CP violation, <p = 0. In this case the relevant mixing angle is 9l — Or, 
see Eq. (4.35). The notation for the various constraints is the same as in Fig. 3. 


with an additional factor N c . As in the case of VLLs, for R w = 1 and large values 
of Q we find two regions of the mixing parameters where /u 77 remains SM-like, while 
can strongly depart from one. (i) For two degenerate masses m\ = m 2 , the 
interference with the SM amplitude is destructive for 77 and constructive for 7 Z, 
as illustrated in Fig. 1. Therefore, there are values of the mixing parameters where 
accidentally /U 77 goes back to the allowed range, while at the same time one can 
even saturate the present upper bound /r 7 z 10. Note that the gluon-gluon channel 
remains nearly SM-like, because its amplitude is not enhanced by the large factor 
Q 2 . (ii) For m\ A m 2 , the amplitude AJJ ^ in Eq. (4.37) can be tuned to zero, while 
at the same time one can have large contributions to /i 7 z together with sufficiently 
small corrections to S and T. E.g. taking N c = 3, Q ~ 8 , mixing parameters ip = 0, 
9l — vr/ 8 , Or ~ vr/10, m\/m 2 — 1.3 and m 2 — 1 TeV, one obtains /r 7 z — 2. For the 
cases R w > 1, we refer to the discussion below Eq. (4.40). 

• The VLQs also contribute to the Higgs production by gluon fusion, with an amplitude 
that can be easily obtained from the 77 one. For the pair of mass eigenstates / 1 , /2 
of charge Q , one has 4® ^ = [3C(R c )/2]/(N c Q 2 )AJJ j 2 , see appendix B. 2 . The 
interference with the SM is constructive in the gluon case, thus enhancing the Higgs 
production. Note that the gluon-gluon channel also receives a non-zero contribution 
from the Q = 0 sector. 

For illustration, we display in Fig. 8 the parameter space for the case R c = 3, R w = 1 
and Y = —1/3 ( Y = —2/3), that corresponds to the pair of VLQs B and Q (T and Q), 
in the limit of no-mixing with the SM quarks. In Fig. 9 we illustrate how the constraints 
change for a larger colour representation, as we replaced R c = 3 by R c = 8 . For definiteness, 
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Figure 9. The same as in Fig. 8, but replacing colour triplets with colour octets. 


we assumed that there is no CP violation, = 0, and that the two mixing mass eigenstates 
fi .2 of charge Q = Y are degenerate in mass, mi = m 2 - In the case of colour octets, 
larger deviations in the Higgs signal strengths are possible, but the various constraints 
are correspondingly stronger. One can reach fi^z — 1.4 for the octets and /r 7 z — 1-2 for 
the triplets. In the case Y = —1/3, /r 77 and are very close to each other, because 
Al h,h/ A Ti,f2 = (9ii + dm) /( Yc w) I s numerically close to one, see Eqs. (4.35) and (4.37), 
and Af M /A^M I s dose 1° one as well. The strongest constraint on the mixing among the 
heavy states comes from the S and T parameters. In the octet case the fit of the main Higgs 
decay channels (see the end of section 5.1 for details) is also a relevant constraint. The 
mass scale mi = m 2 is constrained by the searches of stable coloured particles, discussed 
in section 3. 


• Two VLQs coupled to each other, mixing with the SM fermions. Let us briefly 
discuss the possible interactions between the two VLQs ^2 in Eq. (5.20) and the SM 
quarks. This requires of course R c = 3. A non-zero mixing with the bottom and/or top 
quark can occur if and only if at least one VLQ belongs to the set of seven VLQs in Table 1. 
The complete list is 


singlet + doublet 


T + X T (20 , T + Q (2 1' + b') , 
B + Q{t' + 2b') ,B + Y b {2b') , 


doublet + triplet 


' X T + Z Xt (20 , X r + A q (2 1' + b') , 
< Q + Xq (2 1' + 2b') ,Q + Y q (2 1' + 2b') , 
Y b + Y q {t' + 2b') , Y b + W Yb {2b') , 


triplet + quartet 


Aq + n Xr ( 2 1 ' + 2b') , Xq + n Q {2? + 2b') , 
Yq + n Q {2t' + 2b') , Vq + n YB (2 1' + 2b') , 


(5.21) 


(5.22) 


(5.23) 
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where we indicated in brackets the number of new states mixing with the top and with the 
bottom quark. We also introduced a few new multiplets, with no Yukawa couplings to the 
SM fermions: Zx T ~ (3, 3, 5/3), Wy B ~ (3,3,—4/3), flx T ~ (3,4, 7/6), Qq ~ (3,4,1/6), 
and fiy B ~ (3,4, —5/6). They can be written in components as 


Zx T 



Wy B 



(z\ 


( x\ 


(T\ 

X 

) = 

T 

> ^Yb = 

B 

T 

B 

Y 

\ B ) 


W 


w 


(5.24) 

where the two new exotic states Z and W have charges Q(Z) = 8/3 and Q(W) = —7/3. 
Recasting LHC searches, Ref. [100] puts a lower bound of 940 GeV on the mass of the 
Q = 8/3 state. 

When there is only one t' ( b') state, the mixing in the top (bottom) sector has the 
same pattern as in section 5.1. On the other hand, when there are two t! states, the top 
sector mass matrix takes the form 



(Xt V2 


m 2 \ 

M t = 

m 1 

Mi 





m 2 J 


(5.25) 


where we dropped possible Clebsch-Gordan coefficients. Here ifti (i/j 2 ) is a weak singlet or 
triplet (doublet or quartet), m\ (m 2 ) vanishes unless V’l = T (i/ 2 = Q ), and Ai ,2 vanishes 
if i/j\ 2 is one of the multiplets in Eq. (5.24). The bottom sector mass matrix in presence 
of two b' states has an analog structure. The mixing with the SM is controlled by ratios of 
the type A v/M or m/M. They are mostly constrained by the Zbb couplings and by S and 
T. One typically expects \Xv/M\,\m/M\ < 0.1 — 0.2, in analogy with Figs. 3-6. Possible 
cancellations among the various contributions could relax these constraints. A detailed 
analysis of the whole parameter space of these models is beyond the scope of this paper. 
Some recent study can be found in Ref. [101], that discusses the phenomenology of the 2 1' 
and 2 1' + b' cases. 

Let us focus on the possibility to generate the top mass (and analogously the bottom 
one) through the mixing with the VLQs, in the limit where the SM Yukawa coupling A t 
(A b) vanishes. This is possible whenever the determinant of Ai * in Eq. (5.25) is non-zero 
for At = 0, that is, if and only if the two VLQs both couple directly to the SM. The 
resulting top mass is of order m.t \v(Xv/M) 2 , Xv(m/M) or A v(m/M) 2 . The latter 
possibility is motivated by partial compositeness. In this scenario, the SM fermions do 
not couple directly to the composite Higgs, therefore Xt = Ai = A 2 = 0 in Eq. (5.25). 
Rather, they couple linearly to a composite vector-like fermion with the same quantum 
numbers. This corresponds to the VLQs T and Q for the case of the top quark, leading to 
mt — (my/My) (mQ/M q)Xqtv/ \/2, and analogously B and Q for the case of the bottom. 
The phenomenology of top and bottom partners in composite models, and the associated 
constraints, are analysed e.g. in Refs. [98, 102, 103] (see also Ref. [88] for warped extra 
dimensional models). 
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5.3 Vector-like plus Majorana quarks 


We define a Majorana quark to be a V = 0 fermion multiplet in a non-trivial, real colour 
representation, R c = R c / 1. Such object may couple to the Higgs only in the presence 
of a VLQ in the same colour representation. As a consequence, these new fermions do not 
mix with the SM ones and the lightest mass eigenstate is stable. For the smallest possible 
representation, R c = 8 , the searches for a stable gluino lead to a lower bound ~ 1.3 TeV 
[33]. For bigger representations, N c > 27, one expects an even more stringent limit, given 
the larger production cross-section and a similar hadronisation behaviour. 


• One VLQ plus one Majorana quark. The most general set formed by a Majorana 
quark coupled to a VLQ can be written as 


XR ~ (Rc, Rw , 0) , i/jl, ipR ~ (R c , R w ± 1, -1/2) , R c = R c ^ 1 . 


(5.26) 


If R c is odd, one needs N w / 2 + 4n to avoid the global SU(2) W anomaly. The Lagrangian 
and the structure of the mass matrices are identical to the analogue leptonic case R c = 1, 
see Eqs. (4.46) to (4.52). Here we discuss only the phenomenological differences due to the 
effect of colour. The new states contribute with an additional factor N c to the one-loop 
diagrams for the S and T parameters, as well as for the Higgs signal strengths /r 77 and 
p~,z- In addition, they also contribute to the Higgs production via gluon fusion, with an 
amplitude related to the photon-photon one. For each sector of charge Q that couples to 
the Higgs, one has q = (N C Q 2 ) / [3C (R c ) / 2}A^ l j, q. The gluon-gluon channel receives 
a non-zero contribution even from the Q = 0 sector, that is present when N w is odd. 

Since the new quarks are necessarily heavy, their loop contributions can be estimated 
with good accuracy using the LET approximation, as described in Appendix B.2 and in 
section 4.4. In particular, q is obtained from the amplitude in Eq. (4.51), times a 
factor N c . For the gluon-gluon channel, summing over the different sectors one finds 


^ - k ^C(R c ) 


Xv Xv 
M X M^ ’ 


(5.27) 


where K\ = —4, K 2 = — 2 , A 3 = —8/3, AQ = — 2 , and so on. As the top quark amplitude 
is approximately equal to one, Ajj g gives roughly the ratio between the contribution of new 
fermions and the SM one. Note that either constructive or destructive interference with 
the SM amplitude is possible. The mixing parameters, of the generic form Xv/M, must 
satisfy the constraints from R gg and A 77 , whose allowed ranges are given in Appendix B.3. 
In the minimal case with R c = 8 and R w = 1, only R gg receives a correction, leading to the 
upper bound |Au/M| < 0.17. One expects similar or even stronger bounds from S and T, 
in analogy with the cases of Figs. 8 and 9. Up to possible cancellations, a larger R c leads 
to stronger constraints on the model, and to larger deviations in the Higgs couplings. 


• One VLQ plus two Majorana quarks. Coming to sets formed by two Majorana 
quarks plus one VLQ, the first obvious possibility is to add a second copy of \R to the 
previous case. The additional Majorana multiplet automatically cancels the global SU (2) 
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anomaly, therefore N w is arbitrary. The phenomenology is the generalisation of the one 
discussed above. 

The second and last possibility is provided by the set 


XlR ~ (Rc,Rw,0), fpL, i’R ~ {Rc, Rw + 1, — 1/2), X2R ~ (Rc,Rw + 2,0), R c = R c / 1 , 

(5.28) 

with Ny, necessarily odd if R c is odd. The Lagrangian is the same as in Eq. (4.54), and the 
structure of the mass matrices is also the same as in section 4.4. Let us just present the 
amplitude for Higgs production into gluon-gluon fusion, that is obtained by generalising 
Eq. (5.27), 


A" 

xiX2b 


(Rw) - C(Rc) 



X±v Xiv 


+ K N w +2 


X 2 v X 2 v 

^X2 


(5.29) 


In the minimal case where R c = 8 and R w = 1, the allowed range for R gy leads to a bound 
on the mixing parameters \Xv/M\ < 0.13, where we assumed there is no hierarchy nor 
cancellations among the various mixing parameters. The i? 77 constraint is less restrictive. 


6 Conclusions 

We undertook a systematic analysis of new fermions interacting with the Higgs boson. 
Their properties (gauge charges, masses, Yukawa couplings) are significantly more con¬ 
strained after the measurement of the Higgs mass and couplings at the first run of the 
LHC. It is intriguing to identify the few extensions of the SM that outlived this test. We 
especially aimed at those scenarios that may depart from the decoupling limit, in which 
the new fermions become very heavy and/or their mixing with the SM becomes very small. 

In section 2 we presented the complete classification of sets of n chiral fermions interact¬ 
ing with the Higgs, for n < 4. While the minimal possibilities are well-known, already for 
n = 3 and 4 we singled out several exotic sets of fermions with a peculiar phenomenology. 
They emerge from a non-trivial interplay of several self-consistency conditions: cancellation 
of gauge anomalies, absence of charged massless components, non-zero Yukawa coupling 
to the SM Higgs doublet. In our classification we recovered as a special case the fermion 
content of well-motivated theories beyond the SM, such as the seesaw, supersymmetry, or 
partial compositeness. These cases are situated in a more general context, by considering 
the most general Lagrangian for the new fermions, not restricted by additional theoretical 
considerations. Would the evidence of a new particle emerge from data, one should indeed 
explore the full parameter space, before endorsing a specific model. We also argue that 
larger sets of new fermions, with n > 5, do not allow for qualitatively different phenom¬ 
ena, as all the possible building blocks of a fermion mass matrix already appeared in our 
classification. 

In order to examine the phenomenology of the new fermions, in the appendices A and 
B we derived the general expression of the fermion couplings to the EW gauge bosons and 
to the Higgs boson, for fermions in arbitrary SM representations (R c , R W ,Y). We also 
provided the formalism to define the gauge and Higgs boson couplings to the fermion mass 
eigenstates, after EWSB. Besides these tree-level results, we presented the general one-loop 
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amplitudes for the gauge boson vacuum polarisation, Tly\n, that allows to define the EW 
oblique parameters S and T, and for the Higgs coupling to gauge bosons, hVV', that allows 
to compute the rate for h —>• gg,j'y,'yZ. 

Let us summarise the main results of our phenomenological survey of sections 3 to 5: 

• Several exotic families of chiral fermions, that receive a mass from EWSB only, are 
still marginally compatible with EW precision tests and direct collider bounds. How¬ 
ever, the coloured ones are neatly excluded, as they would greatly enhance the hgg 
coupling. On the other hand, a colourless family formed by two weak doublets and 
four singlets is still compatible with the measured /177 coupling. 

• The mixing of two or more sterile neutrinos with the SM leptons can have observable 
effects, despite the smallness of the neutrino masses. If the sterile neutrinos are 
lighter than the EW scale, they may modify significantly the Higgs invisible width; 
if heavier, they can appreciably contribute to the S and T parameters. 

• In general, a heavy charged lepton t’ cannot mix significantly with the r because 
of the ZTT-coupling constraint. Nonetheless, in a few special regions of parameters 
interesting phenomena are possible: (i) If m T i < mh, the decay rate for h —> t't can 
be significant despite the small mixing, (ii) When both tr and tr mix with heavy 
leptons, it is possible to generate m T entirely through the small mixing permitted 
by the Z couplings, as long as the two heavy leptons are connected by a Yukawa 
coupling A > 3. (iii) If the new lepton sector is arranged to have an approximate 
custodial symmetry, the Z couplings are protected. In this case a large r — t' mixing 
is allowed, and it may strongly suppress h —> tt. 

• There are two extended classes of new fermions that can couple to the Higgs dou¬ 
blet without involving the SM fermions: either a pair of vector-like fermions, whose 
components can have an arbitrary charge Q, or a pair formed by a vector-like and 
a Majorana fermion, whose components have (demi-)integer Q. These fermion pairs 
were not studied in full generality in the previous literature, and they can produce 
large observable effects, even when the mixing with the SM fermions is zero. By vary¬ 
ing their mass matrix parameters, one can typically scan over the full allowed range 
for the signal strength // 77 , while remaining in agreement with direct collider searches 
and EW precision tests. In most cases receives a correction comparable to /z 77 , 
but when the latter is accidentally close to one, it is possible to have 5g^z <5/Wr 
We will discuss this point in detail below. 

• The mixing 6 of a heavy quark t' or b' with its SM partner is constrained by the 
EW precision tests. Nonetheless, the b — b' mixing may significantly suppress the hbb 
coupling, leading to corrections as large as d/r 77 ~ dg^z — 0.6. We also notice two 
remarkable circumstances that allow for a large mixing: (i) The corrections to Zbb 
from t — t! mixing are loop-suppressed and may be also suppressed by (mt/m t ') 2 - The 
T parameter receives opposite sign corrections that cancel each other for a specific 
value of sin (9 x mg. Both conditions can be realised in the case of the VLQ doublet 
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(A, T), allowing for a large sin# < 0.5. (ii) The Lagrangian of the new quarks can 
preserve a custodial symmetry, that suppresses the corrections to T as well as to the 
Zbb couplings. In this custodial subspace of parameters, the upper bound on the 
mixing relaxes, the exact value depending on the model: for the VLQ doublet Q 
coupled to t.R and we find sin # < 0.15, while for the two doublets Q + Yb coupled 
to bji (Xt + Q coupled to f#) one can reach sin# < 0.45 (sin# < 0 . 8 ). 

In the course of our analysis, we paid special attention to the relative contribution 
of the new fermions to h —> 77 and h —> 7 Z, as the former rate is already constrained 
to be close to the SM prediction, while the latter could still depart strongly from its SM 
value. It is commonly believed that new physics cannot provide a large correction to the 
7 Z channel without affecting 77 as well. Indeed, let us consider the effective Lagrangian 
before EWSB, that corresponds to the limit where the new fermions are heavier than the 
EW scale. There are several dimension-six operators involving the Higgs doublet H and 
the field-strengths BWf w . listed e.g. in Ref. [104], The operators contributing to / 17 Z 
can be generated, at one loop, only by two fermion multiplets coupled to H. At least 
one of these fermions has non-zero hypercharge, thus it necessarily induces the operator 
W as well, that contributes to h. 77 . One can rephrase the same argument in 
terms of the effective Lagrangian for the hVV' couplings after EWSB, that is displayed in 
Eq. (B. 8 ). The coefficients of the dimension-five operators, generated at one loop by the 
fermion mass eigenstate fi, are given in Eq. (B.22) for /177 and in Eq. (B.34) for L 7 Z. 
The fermion fi cannot contribute to the 7 Z channel only, simply because one needs a 
charge Qi ^ 0 and a non-zero coupling y* (or yj) to the Higgs, therefore the 77 channel 
receives a contribution too. This argument, however, has some loopholes: first, the sum 
over all fermion mass eigenstates can lead to a cancellation in the signal strength /j 77 and 
not in y.^Zi as the summands in the two channels differ by a factor ~ gj' /Qg, second, / 17 Z 
receives an additional contribution from loops involving two fermion mass eigenstates, with 
off-diagonal couplings to both h and Z, see Eqs. (B.31) and (B.32). 

As a matter of fact, in our survey of fermionic extensions of the SM, we encountered 
a few scenarios where <5// 7 z 3> hy 77 : 

(i) One can exploit the order one differences between the Z and 7 couplings and loop 

functions. For example, in the case of t — t' mixing, by 77 is proportional to A 1 / / 2 (r/) — 
A 1/2 (77), that is very small as both form factors are close to the asymptotic value 
Ai/ 2 (0). On the contrary, the correction to y^z is controlled by Ay) — 

y^A 1 / 2 (rt, A r ), that is in general of order one. Also, off-diagonal loops provide an 
additional contribution of the same order. Unfortunately, the absolute size of the 
correction is too small to be observed, as the mixing between the SM and new fermions 
is subject to the EW precision constraints. We find at best hy 77 <C Sg^z — 0.2, hardly 
visible even with 3000 fb -1 at 14 TeV, see Table 5. 

(ii) A much larger 5y^z is possible when new fermions couple to each other through the 
Higgs. Each sector of heavy states with given charges N c and Q gives a contribution 
to the /177 amplitude proportional to yiv/rrii. We found that the structure of the 
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fermion mass matrix allows this sum to vanish, see the discussion below Eq. (4.39). 
At the same time, the 7 Z amplitude is proportional to JV g(yiv/rrii. One can obtain 
e.g. n 1 z — 2, by means of a pair of states with N c Qs ^ ~ 3, where s^ is the relevant 
mixing parameter. Alternatively, the same effect is produced by several states with 
smaller charges. The required set of parameters can be in agreement with S and T 
as well. This opens a discovery opportunity for the second run of the LHC. 

(iii) There is a second possibility to achieve a large y 1 z- The signal strength fi 77 may be 
accidentally close to the SM, because the amplitude generated by two new fermion 
multiplets coupled to the Higgs has sign opposite to the SM one, and for Aj 1 ~ 
— 2 A'gM one recovers // 77 ~ 1. One needs either small weak multiplets with large 
charges, N c Q 2 s ^ ~ 5, or larger multiplets with smaller charges. In this region of 
parameters the S and T constraint can be satisfied and, moreover, one generically 
expects fi^z much larger than one, because for large values of Q the amplitude Aj Z 
interferes constructively with the SM. We find that one can almost saturate the 
present experimental bound y^z 1$ 10, therefore coming LHC data will be able to 
quickly probe this scenario. 

In the cases listed above, a mild tuning of the parameters is sufficient to comply with the 
presently allowed range for /i 77 , shown in Table 4. In the future, the room for a large // 7 z 
will progressively shrink. 

The second run of the LHC, that recently started data taking at 13 TeV, will close in 
on most of the scenarios we have been considering. The allowed regions of parameters at 
low masses will be covered by direct searches for new fermionic resonances. The islands 
that survive at large mixing between the SM and new fermions will be probed by the 
increasing precision in the Higgs coupling measurements, even though there are models 
where one needs to wait for a high accuracy. In the absence of a signal, we shall be 
virtually cornered to the region of very heavy masses and/or very small mixing. Even 
when the new fermions are too heavy to be directly produced and mix negligibly with the 
SM, their Yukawa couplings to the Higgs can be effectively constrained by the radiative 
Higgs couplings. 
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A Electroweak precision tests in presence of new fermions 

In this appendix we provide general formulas for the EW gauge boson couplings to fermions, 
as well as for their vacuum polarisation amplitudes. This allows to define and compute the 
oblique parameters S and T [105-108]. We discuss the experimental constraints on these 
parameters, as well as on the Z couplings to light SM fermions, such as Zbb. 


A.l Electroweak gauge boson couplings 

The couplings of the EW gauge bosons to a chiral fermion multiplet in a given representa¬ 
tion (. R W ,Y) of the EW gauge group are determined by the covariant derivative 

D, = d„ - igT a W; - ig'YB „ = d„- iJL(T+W+ + T~W~) - igT 3 W 3 - ig'YB „ , (A.l) 

where T a are the SU(2) W generators for R w , T ± = T 1 ± iT 2 and = (W* zp iW 2 )/\/2. 
In full generality, the resulting non-vanishing couplings are 

C W± q^,U = = %ViU+_ 1) - W ± , (A 2) 

c w3 (fm',fm) = gT^5 m ' m , = g'Y5 m ' m , 


where j = {N w — l)/2 is the weak isospin, and the N w components of R w are labelled by 
m, m' = —j, —j + 1,..., j — 1, j, and have electric charge Q m = Tf n +Y. It is straightforward 
to derive from Eq. (A.2) the couplings of = c w W 3 — s w B M and = s w W 3 + c w B 

C (/m'j/m) (Zrri j C ( fm'ifm .) &Qm$m'm • (A.3) 

C-w 

After EWSB, for each value of the charge Q , the fermion mass term can be written as 
fLa(-M- q) aft fRj3 1 where a, (3 = 1,..., hq run over the uq fermions of charge Q (in a given 
colour representation). In general the mass matrix is not diagonal and the mixing can be 
described by f La = (i UfydfLi and f Ra = ( UQ) ai f Ri , where ./) are the mass eigenstates. 
Therefore, the couplings of the gauge bosons to the mass eigenstates are 


(<*)« = ^ = W ± ,W 3 ,B , (A.4) 

where the (c\ R ) a /3 are given in Eq. (A.2), and Q = Q' ± 1 for V = W ± , Q = Q' for 
V = W 3 ,B. 

The mixing cancels out in the photon couplings, because U(l) em is unbroken, and 
one finds immediately (c A R )ij = eQ5ij. The Z-boson couplings to the mass eigenstates, 
instead, do depend on the mixing. Using the parametrisation 


= f z f E («.A +u = fi->“ (»«- 4U /. 

(-"w . . t- w ■ ■ 

h3 h3 

one finds 

alf = Tl8 o3 {U^ R )Y(U^ R )s 3 - slQSij , g v > A = 


j > 


g L ± g R 


(A.5) 


(A.6) 
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The matrices g L,R (g v,A ) are hermitian, with possibly non-vanishing off-diagonal entries. 
Note that the mixing of fermions with equal EW charges does not affect the couplings to 
the neutral gauge bosons: if T 3 (or, equivalently, Y a ) is the same for all a, then one can use 
{UQ R )* ai {UQ R ) a j = Sij. and the couplings to Z (as well as to W 3 and B ) reduce to their 
unmixed values. Thus, the neutral-current couplings of SM fermions receive a correction, 
only when they mix with new fermions with a different value of T 3 (of Y). For example, 
the mixing of two left-handed fermions fi, a and fn, of charge Q amounts to 



9 L = 


T 3 - 

1 a 


4 Q 


n 3 - 


s wQ, 


+ ( T a - T b) 


— S 


SC 


In this paper we assume that new fermions mix with the third SM family only. Indeed, 
flavour-changing neutral currents among the different SM families are strongly constrained 
experimentally. 


A.2 Constraints from S and T 

The vacuum polarisation amplitudes for the EW gauge bosons, defined by the effective 
momentum space Lagrangian 

£n = -W+n% w (p)W- - l -B,n^(p)B u - W 3 U^ 0 (p)B u - -w;;u^(p)w 3 , (A.7) 

can be decomposed into transverse and longitudinal parts, 
u vv'(p) = n y-V' (p)(p 2 fj ll,y -P ,J 'P U ) + n vv' (p)jAA = TI VV > (p 2 )g lw + {p ll p v - terms) . (A.8) 

As in the experiments the mass of the external fermions is much smaller than the EW 
scale, rrij <C m 2 z , one can drop the p^p^-terms and expand in p 2 , 


n V v{p 2 ) = n VV r(o)+ P 2 u' vvr (o) + o(p 4 ). (A.9) 

The lowest terms in this expansion are sufficient to describe accurately the effect of heavy 
new physics: when new particles at scale uif contribute to the vacuum polarisation am¬ 
plitudes, the higher order corrections are suppressed by powers of m 2 z /m 2 F . Taking into 
account that three coefficients can be traded for the experimental values of a, s w and 
mz, and two others are determined by the Ward identities for the photon, one finds that 
two parameters are sufficient to characterise the effect of new physics at leading order in 
m 2 z lm 2 F [107, 108]. The combination that describes the custodial symmetry breaking at 
leading order is given by 

t = [(n 33 (o) - nf 3 M (o)) - (iw(o) - n^(o))] . (a.io) 

ac^m z 

The combination that breaks the weak isospin at leading order, but respects the custodial 
symmetry, is given by 

5 [(n 30 (m|) - (n 30 (o)) - (nf 0 M (m|) - nf 0 M (o))] ~ ^^[n' 30 (o) - n' 3 s 0 M (o)]. 

arriy a 

(A.ll) 
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Figure 10. The 68 % (red), 95 % (orange) and 99 % (yellow) C.L. ellipses in the S — T plane, from 
the fit of Ref. [109], with the other EW parameter U left free. The black dot indicates the best fit, 
while the star at S = T = 0 is the SM point, with mt t ref = 173 GeV and mh,ref = 125 GeV. 


fi 



Figure 11. Contribution of two fermion mass eigenstates /i and fa to the vacuum polar¬ 
isation amplitude for the gauge bosons V^ and V ' v . The relevant couplings are defined by 
Cf fv - vjtrMhPL + (cDijPRifj. 


The approximation in terms of amplitude derivatives evaluated at p 2 = 0 is appropriate 
only for new physics much heavier than mz that does not mix with light SM particles; 
in the general case one should keep the definition of S in terms of amplitude differences, 
to avoid unphysical singularities that may appear in the derivative. The subtracted SM 
contribution is evaluated at a reference point for the SM parameters. Following Ref. [109], 
if one takes mt, re f = 173 GeV and mh t ref = 125 GeV, the present experimental allowed 
ranges are given by 

S = 0.05 ±0.11, T = 0.09 ±0.13, (A.12 ) 

with a correlation coefficient ~ 0.9. In Fig. 10 we display the allowed region in the S — T 
plane, that we adopt in the rest of the paper to constrain the parameter space of each 
model. 

In order to estimate the contributions to S and T, in any theory where the new physics 
is weakly coupled, one should just compute the one-loop diagram contributing to the EW 
gauge boson vacuum polarisation amplitudes, shown in Fig. 11. The most general couplings 
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of the EW gauge bosons to the fermion mass eigenstates are defined in Eq. (A. 4). The 
functions Ylyy/(p 2 ) defined in Eq. (A. 8 ) will receive different contributions from the left- 
and right-handed couplings, 

n w = eld Ull + dd H.R-R + ddi H lr + did n rl ■ (A.13) 


Performing the computation with dimensional regularisation, one finds the following gen¬ 
eral result: 


n ll,rr(p 2 ) = 


N r . 
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77117712 2 


47T 2 /1 
- 7 

77117712 O 


R(p 2 ) (A-14) 
(A.15) 


where N c is the dimension of the SU(3) C representation of the fermions / 12 , my 2 are the 
masses of / 17 , e and p are defined by d A k = p 4 ~ n d n k and 1/e = 2/(4 — n) — 7 — ln 7 r with 
7 ~ 0.5772, and finally 


2 s _ l/A 771 2 + m 2 — a/A — (p 2 + ie) 


R(p ) = ~y log ■ 


2771.17712 


A = (ml — ml) 2 — 2p 2 (m 2 + ml) + p A . 

(A.16) 

Note that R(p 2 ) is invariant for —7 — \/A and it has a non-vanishing imaginary part 
for p 2 > (? 7 ii + m 2 ) 2 • When evaluating S in Eq. (A.11), one must include only the real 
part of R(m 2 z )i while the imaginary part contributes to the decay width of the EW gauge 
bosons, that may be also modified with respect to its SM value. Note also that, when /1 2 
are Majorana fermions, one must include in the amplitude an additional symmetry factor 
1 / 2 . 

Let us provide for convenience the form of II yy(p 2 ) in some relevant limits. For 
777.1 = in 2 = 0 one has 


tt r 2\ 21V C 
73ll,rr(p ) = 


- +log 
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4ttV 7 

— (p 2 + ie) 6 


3 ( 477 ) 2^^6 + 1 ° g -(p 2 + ie) + 6j ’ U lr,rl(p 2 ) ( A1 ‘) 

This limit is relevant for loops involving the light SM fermions (all but the top quark), whose 
mass can be neglected. In order to compute T and S (in the derivative approximation), 
it is sufficient to compute the first and second term in the p 2 -expansion of II yy(p 2 ), 
respectively: 
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(A.18) 
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For mi = 0 and 777-2 = m (relevant e.g. for the bottom-top quark loop) one obtains 
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= 3(4^ (? + ln -HT - 6, 

(A.20) 

and IIi7j j ^i(0) = U' LR RL (0) = 0. For m± = 777-2 = m^0 (loops involving a unique fermion 
mass eigenstate) one reduces to 
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(A.21) 
(A.22) 


As an illustrative example, consider the case of a fermion “family” with no mixing 
with the SM fermions, formed by one weak doublet Ql = (Tr, Bj ) ~ (R C ,2,Y) and two 
singlets Tr ~ (. R c , 1, y + 1/2) and Br ~ (. R c , 1 ,Y — 1/2), with dim(i? c ) = N c . After EWSB 
they combine into two mass eigenstates T and B with masses 7775 and 777 , 5 ; their non-zero 
couplings to EW gauge bosons are obtained from Eq. (A.2), 


^ ± (T,B) = S-, cT(T) = f, cjT(B) = -f, 

V f l\ ( l\ (A. 23 ) 

c?(T) = c?(B) = g'Y , c B R (T)= g '[Y + - J , c£(S) = </(E - - J . 

To compute the correction to S, one should evaluate Eqs. (A.14-A.15) for p 2 = 777 ^ and 
777i = 7772 = 777 j 1 5 , while for p 2 = 0 one can use directly Eq. (A.21) with m = mR y B- Adding 
the various contributions to II 30 as shown in Eq. (A.13), and replacing into Eq. (A. 11), the 
result is 


(A.24) 

If one adopted the approximate expression for S in terms of derivatives, given by the right- 
hand side of Eq. (A.11), then using Eq. (A.22) one finds only the first term in the squared 
bracket of Eq. (A.24), which is accurate for 777^5 3> mz- To compute the correction to T, 
one should use Eq. (A.21) for the T and B loops that contribute to II 33 , and Eq. (A.18) 
for the T/B loop that contributes to Tlww- Replacing into Eq. (A. 10) one obtains 
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(A.25) 


Particularising these results to the case of the SM top and bottom quarks ( N c = 3, 
Y = 1/6), and neglecting the uncertainty on 7777 , as well as ( 777 ^/? 77 ^)-corrections, we can 
immediately extract the well-known dependence of S and T on the value of the top quark 
mass, 
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A.3 Constraints from Zff 

The Z-boson couplings to the SM fermions are precisely measured. We discuss only those 
with the third family, since in this paper we assume that the mixing of the new fermions 
with the light families is negligible. The deviations with respect to the SM can be expressed 
in terms of the Z partial decay width into any given final state //, 

R(Z -> //) = T{Z ~* fl) = 1 + SR(Z -> //) , f = u r ,T,b. (A.27) 

V SM {Z ^ //) 

L R 

For each fermion / there are two independent couplings g^ as shown in Eq. (A.6), that 
can be separately constrained if the angular distribution of the fermions is measured. 

Beginning from leptons, the Z invisible width and its width into taus are determined 
at the per mil level [3], T(Z —» inv ) = 499.0 ± 1.5 MeV and T(Z —> r + r”) = 84.08 ± 0.22 
MeV. Given this precision and the relatively good agreement between the central values 
and the SM predictions, we constrain the mixing with the new leptons by imposing a rough 
3cr upper bound, 

\8R(Z->inv)\ < 9- 1(T 3 , \5R{Z t + t~) | < 8 • 1(T 3 . (A.28) 


Coming to the Z coupling to bottom quarks, a more detailed discussion is worth, to 
fairly gauge the resulting constraint on the new fermions. The Zbb Lagrangian can be 
written as 


P-Zbb ~ ~Z^b 


{sbb,SM + ^9bb) A ( 9bb,SM A ^9, 


Pr 


(A.29) 


where the SM couplings at tree-level are given by gj^ SM = —l/2 + s^/3 and g *| SM = s^/3. 
Deviations at tree-level occur when the bottom quark mixes with a new fermion with a 
different value of T 3 (and Y). The present experimentally allowed range is given by [110] 


5g A = 0.0016 ± 0.0015 , Sg$ = 0.019 ± 0.007 , (A.30) 

with a correlation coefficient ~ 0.8. In Fig. 12 we display the allowed region in the 
plane. Note that the left- and right-handed couplings are determined with per mil and per 
cent precision, respectively, and the best fit region is incompatible with the SM at about 
99% C.L.. In some analyses, a slightly better agreement is obtained, at about 95%, due 
to different details in the global electroweak fit, see e.g. Ref. [111]. The discrepancy with 
the SM comes mostly from the measurement of the forward-backward asymmetry A b FB , 
that may be just an upward statistical fluctuation, an unidentified systematic error, or 
alternatively an indication for a significant new-physics contribution. In this paper, in 
order to constrain the new fermions that modify the Z&5-couplings, we will conservatively 
enlarge the 99% C.L. region, by allowing it to shift towards the SM point, till the latter 
touches the 68 % ellipse, as illustrated in Fig. 12. Such a shift roughly corresponds to 
introduce a systematic error in the measurement of A b FB . 

The mixing with new fermions can be such that uo tree-level deviations occur in the 
Zh 6 -couplings, but they do occur in the Wtb coupling and/or in the Ztt couplings. These 
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Figure 12. The 68 % (blue), 95 % (magenta) and 99 % (cyan) C.L. ellipses in the 8g *| — 8g F 
plane, extracted from Ref. [110]. The black dot indicates the best fit, while the star at the origin 
represents the SM. In our analysis, we allow for a larger parameter space, delimited by the dotted 
line, that is obtained by shifting the 99% ellipse towards the origin, till the SM point enters the 
68% region. 


L R 

deviations may affect significantly , because the contribution of one-loop diagrams 
involving the top quark and the TY-boson is larger than the experimental uncertainty. Also 
new fermions may correct significantly Zbb at the one-loop level, if they are not much 
heavier than the top. The detailed structure of the one-loop corrections to Zbb can be 
found e.g. in Ref. [112], that we employ for our analysis of models with modified top 
couplings. For example, the correction to the left-handed coupling, from the top loops 
only, can be written as 

= 5gttf L {m W /mt) + 8g R f R (m w /m t ) + Sg^f W (m w /m t ) , (A.31) 

where bg^ R are defined in analogy to Eq. (A.29) with t o b, the Wt^L coupling is given 
by (g/V 2)(1 + bg^), and the functions / l,r,w can be extracted from Ref. [112], 

B Higgs boson couplings in presence of new fermions 

In this appendix we present a general parametrisation for the Yukawa couplings among 
the SM Higgs doublet and two arbitrary fermion multiplets. We then analyse the resulting 
modifications in the Higgs couplings to the SM particles, at leading order: corrections 
at tree-level to the Higgs-fermions couplings, and one-loop corrections to the Higgs-gauge 
bosons couplings. Finally, we briefly review the present experimental constraints on these 
couplings. 

B.l Tree-level Higgs couplings 

The SM Yukawa couplings are given by 

- £y M = VuCIlurH + y d qLd R H + y e l L e R H + h.c. , (B.l) 
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where qL = ( ULdi) T , H = (H + H°) T and H = ia-zH* are SU(2) W doublets, while ur, cIr 
and eR are singlets, and flavour indexes are understood. 

In full generality, the Higgs doublet may have a non-zero Yukawa interaction with any 
pair of chiral fermions that transform under SU (3) c x SU(2) w x U (l)y as ipL ~ (-R c , Rw, Y) 
and ip% d ~ (R c , R w ~ 1 , Y ± |): 


- C Y = Y 

4>l 


Y (mrh) +X] 

.n k 


+ h.c. . 


(B.2) 


Here the parentheses stand for the appropriate contraction of the SU(2) w indexes. Let us 
denote the components of ipL by the index m = j, j— 1,..., — j+1, — j, where j = (N w —l)/2. 
Then, the multiplet with components (V’L)m = transforms in the conjugate 

representation R* w . It is possible to define a multiplet V’L that properly transforms in the 
representation R w , by using the SU(2) w conjugation matrix R, (ipL ) m = R mn ('^L)n = 
(— iy~ m (il>L)-m- The 2 j components of pair with the upper (lower) 2 j components 
of ■) Pl to form the upper (lower) component of a weak doublet D. The corresponding 
Clebsch-Gordan coefficients are given by 



1 1 

-;m,-m± - 



±(-l> 


3-m 


I j ±m 

j(2j + 1) 


(B.3) 


Contracting D with H into a weak singlet, D a (ia 2 ) a bHbi one finds 
3 


Y 


m=-j-\-l 


I j + m 

J{ 2j+l) l 


• (B.4) 


The same expression holds for ip^ ip 1 ^ and H <->• H as well. Of course, all the results 
above also apply when one makes everywhere the replacement L -H- R. 

The relative size of the Clebsch-Gordan coefficients has important phenomenological 
consequences, e.g. the different components of the fermion multiplets acquire a different 
mass after EWSB. The overall normalisation of the SU(2) w contraction is also important, 
to establish the perturbative range for a Yukawa coupling y: for instance, the contribution 
of y to the Higgs wavefunction renormalisation at one-loop goes as y 2 /(167r 2 ) times the 
sum of the Clebsch-Gordan coefficients squared, taken over all possible isospin components 
in the loop. Adopting the above conventions, such a sum is normalised to one, and we can 
easily define the region where perturbation theory can be trusted, by requiring y/ (47r) <C 1. 
However, we keep the conventional normalisation for the doublet-doublet contraction into 
a singlet, with no overall factor 1 /\/2, that strictly-speaking should be included: the issue 
of SU(2) W normalisation is more relevant for large weak multiplets. 

We note that the perturbative upper bound on a Yukawa coupling y depends on the 
process under consideration. Schematically, the next-to-leading order amplitude is given 
by the leading order one times a factor y n g m F c / where g stands for other couplings 

such as gauge couplings, and F c is the colour factor, with typical values F c = 1, N c ,C(R c ). 
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In the example of the Higgs wavefunction normalisation adopted above, one has n = 2, 
m = 0 and F c = N c , therefore we could have taken into account the SU( 3) c contraction 
by adding a factor 1 /\fN~ c on the right-hand side of Eq. (B.4), or alternatively requiring 
yy r Nj (4vr) <C 1. However, the one-loop amplitudes relevant in our analysis (EW precision 
tests, Higgs couplings to fermions and gauge bosons, etc.) behave differently from each 
other, and the perturbativity criterion varies correspondingly. In some cases the colour 
enhancement is absent, or compensated by small gauge couplings, or by a small mixing 
between new and SM fermions. Therefore, we find more conservative to stick to the bound 
y -C 47 t. 

At EWSB, H° can be replaced by ( v + h)/\/2 , to obtain the couplings of the phys¬ 
ical Higgs boson h to the fermions in the interaction basis. All fermions with equal 
charge and in the same colour representation may mix, and their mass matrix A4 = 
UL,diag(m\,... ,m n )U\, may include both v-independent vector-like mass terms, and the 
EWSB contributions ~ yv. Thus, one can derive the /i-couplings to the fermion mass 
eigenstates as follows: 
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(B.5) 


where y and y are hermitian matrices defined by 
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A +At 
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A-At 
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(B.6) 


In the CP-conserving case A is real, therefore y = y T is real and y = —y T is imaginary. In 
the case of purely chiral masses (e.g. in the SM), one has M(v) oc v, therefore dM.{v)/dv = 
M/v, y = A = diag(mi,... ,m n )/v and y = 0. On the other hand, in presence of both 
chiral and vector-like masses, the Higgs boson can have both scalar and pseudo-scalar, CP- 
even and CP-odd, diagonal and off-diagonal couplings to the fermions mass eigenstates, 
and its couplings are not proportional to the fermion masses. A simplification occurs 
in those SM extensions such that (dM/dv) a p = M. a pcp, that is, each row of the mass 
matrix has the same dependence on v. In this case t/j) and U k cancel out in A and one 
finds X jk = m j J2i3Cp(Uf t )pj(U R )p k . Similarly, when ( dM/dv) a p = c a M a / 3 , one finds 
A jk Cn(f^r)ni ( Ulj)nk■ 

The mixing with new fermions modifies the Higgs boson decay width into SM fermions 
at the tree-level. In addition, there may be new Higgs decay channels, with one or more 
new fermions in the final state, as long as they are lighter than h. In full generality, the 
Higgs decay width into two fermions at leading order is given by 
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9{m h - rrij - m k ) , 


(B.7) 
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where /3± k = [1 — (rrij ± mk) 2 /m 2 ] 1 / 2 , Ajk = 2 if the final state particles are identical 
Majorana fermions (j = k and fj = /?), and A jk = 1 otherwise. 

Besides the Yukawa couplings to fermions, the Higgs boson has tree-level couplings 
to WW, ZZ and to itself. In the presence of new fermions, all the tree-level couplings 
may receive corrections at the one-loop level. Even though these corrections may become 
relevant in view of precision measurements of the Higgs couplings, they represent in general 
a sub-leading effect and we will not discuss them further. In the following subsection we 
will focus instead on a more sensitive probe of new physics: those Higgs couplings that are 
absent at tree-level in the SM. 

B.2 Loop-induced Higgs couplings 

At one-loop new couplings are induced between the Higgs boson and the SM particles, 
that are absent at tree-level. In particular, U(l) em and SU( 3) c gauge invariance prevents 
renormalizable couplings to photons and gluons, therefore the tree-level amplitudes for 
h —> gg , h -> 77 and h —>• Z 7 are zero, and in addition the one-loop amplitudes for 
these processes are free from divergences. The effective Higgs boson couplings generated 
at one-loop can be described in full generality by two dimension-five operators, 

C-hvv = ^c/iVv r 'H AU ,H ,AU ' + -ChW'VnvVpcr e ^ wpCT ' S J h > (B.8) 

where W = gg, 77,7Z and V /w = G a ^ u ,F^ v , Z gv are the field strength tensors for the 
gluon, the photon and the Z, respectively. The CP-even (odd) coefficients Chvv {phW’) 
have mass dimension minus one, and may receive contributions from loops of both SM and 
new particles. One should be careful to avoid double-counting: if one considers C-hW as 
part of the effective Lagrangian valid below the EWSB scale v, Chvv and ChVV' should 
include only the contributions of particles heavier than v. 

In momentum space, the hVV' couplings are given by 

£h.W'(p,p') = 2 [chw{p ,fl p v ~ P ■ p'g pu ) + ChW'^^PpPa] h(p + p')Vix(p)V'v(p') • (B.9) 
The decay width of the Higgs boson into two vector bosons is then given by 

+ Oi iW'Pvv'Pvv i 

(B.10) 

where N c = 8 (1) for gluons (for 7 and Z ), A yv' = 2 for V = V' and Ayy = 1 for V ^ V ', 
and finally Py V , = [1 — (my ± my/) 2 /m|] 1//2 . In the following, we will match the explicit 
loop computation with the effective coefficients cy t yy an d Chvv' ■ 

B.2.1 Higgs coupling to two gluons 

The Higgs-gluon-gluon coupling Ch gg ( Ch gg ) receives a one-loop contribution from each 
coloured fermion with a non-zero CP-conserving (CP-violating) Yukawa coupling to the 
Higgs boson. The gluon coupling to fermions is determined by SU( 3) c gauge invariance, 

£gff = Jib^A“(T“) bc f ic , (B.ll) 
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Figure 13. Fermionic triangle loop contributing to the coupling of the Higgs boson to two gauge 
bosons. A crossed diagram with the two gauge boson insertions interchanged has to be added too, 
and one must sum over all possible sets (/,, fj,fk ) of fermion mass eigenstates. For the gluons one 
has G a ik = g s SikTf, for the photon G] k = eS ik Qi , and for the Z boson Gf k = g(gY k - 75 gtk)/ c w 


where T“ are the SU( 3) c generators in the representation R c i of the fermion f t . As the 
SU(3) C symmetry is unbroken, there are no ‘off-diagonal’ gluon couplings to two different 
mass eigenstates. Recall that the Higgs is mostly produced via gluon fusion, with a partonic 
cross-section that can be expressed as a function of the partial decay width, 
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8 m h 

In the SM, the contribution of the quarks triangle loop to h —>• gg reads 
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(B.12) 
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where r q = m^/(Am q ) and the form factor is given by 
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(B.14) 

As illustrated in Fig. 14, the top quark gives the dominant contribution, because r* <C 1 
and Th e,... 1- In a generic extension of the SM, the fermions will couple to the Higgs as 

in Eq. (B.5), but only the diagonal, real couplings y k = yu and y k = fja are relevant for 
h —>• gg. One obtains 
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where the CP-even and CP-odd amplitudes are 
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Figure 14. Real (blue lines) and imaginary (red lines) parts of the form factors Ai/ 2 (t) (solid lines) 
and Ai/ 2 {t) (dashed lines). The horizontal lines correspond to the asymptotic values 2, 4/3 and 0. 
The vertical lines correspond to the reference values T t — 0.13, r = 1, Tf, ~ 230 and t t ~ 1300. 


with the Dynkin index C(P c ,) defined below Eq. (2.1). 

By matching with Eq. (B.10), one finds that the contribution of a fermion loop to the 
effective / 155 -couplings is 
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(B.17) 


where c hgg = \ J2i c h gg I and Q hgg = \ Yli&hgg I- In t]i e heavy fermion limit, 2 m, > rn h , one 
can use ^ 1 / 2 ( 0 ) = 4/3 and A 1 / 2 (0) = 2. We note that, in the literature, a factor 1/2 is 
sometimes missing in the expression for c^gg- 

One can use the Low Energy Theorem (LET) [113, 114] (see also [115-117]) to evaluate 
the effective hVV' couplings induced by states much heavier than the EW scale. For a given 
sector of mixing states in the representation ( R C ,Q ) of SU(3) C x £7(l) em , the low energy 
result is a function of their mass matrix Ai only. For the CP-conserving and CP-violating 
[118] gluon-gluon case one finds, respectively, 


LET 
u hgg 


T 


det 


c h gg = ^C(P c )^ arg[det(M)] . (B.18) 


This is very useful in the case of a large, complicated mass matrix A4, because this ex¬ 
pression is much easier to evaluate, with respect to an explicit computation of the mass 
eigenvalues rrq and of the mass eigenstate couplings y, and y,;. Note, however, that this 
approximation requires all the mass eigenstates in a given sector to be heavy, 2 m, 2 > rrih- 
It is easy to check the consistency of Eq. (B.17) and Eq. (B.18) for one heavy chiral fermion 
(e.g. the SM top quark), as = mi = y,u and ^4i/ 2 (p) — 4/3. 
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B.2.2 Higgs coupling to two photons 

The fermions charged under U(l) em contribute to the Higgs-photon-photon couplings c/i 77 
and c hry at one loop. In the SM, there is also the contribution from TT-boson loops, that 
we include using the SM tree-level couplings of the W to the Higgs and to the photon. The 
SM decay width is given by [113, 114, 119] 


TsmO -1 77) = 


2 S 

a mz 


A\{jw)+ ^2 N C iQjAi/ 2 {Ti) 

heSM 


(B.19) 


2567r 3 u 2 

where T\y = m|/(4m^), r* = m|/(4m|), and the form factor for the IT-loops reads 

_ _2^> + 3 t + 3(2t-1)/ W 

T z 

The W contribution is dominant, A\ {t\y) ~ —8.36, and it interferes destructively with the 
top-quark loop, N ct Q^Ai/ 2 {Tt) — 1.83. 

In presence of extra fermions, there are new contributions to h —> 77 that depend, as 
in the case of h — > gg, on the Yukawa couplings y* and y t , 
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In terms of the coefficients of the effective Lagrangian Eq. (B. 8 ), the contribution of each 
fermion is given by 
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(B.23) 


The LET approximation, for a set of heavy fermions in the representation (R c , Q) with a 
mass matrix A4, is given by 


LET _ _?L n 2 N A l n 
Ch77 “ 12 tt Q c du 


det 


Cfe 77 ’ = 2zQ' z N c -^- arg [det (A4)] , (B.24) 


47T 


in analogy with Eq. (B.18). 


B.2.3 Higgs coupling to a Z boson and a photon 


The last loop-induced coupling to be considered is hZ 7 . It is generated by IT-boson loops, 
that we take to be SM-like, as well as by fermionic triangle loops. The Z-boson couplings 
to fermion mass eigenstates are defined in Eq. (A. 6 ). 

The Higgs decay width into a photon and a Z in the SM is given by [120, 121] 


r SM^h —> 7 Z) 


2 2 S / 
Qfl c i v m h ( 1 

5127t 4 u 2 y 



a.i(tw,\w)+ 

Si&SM 


NciQid, 


V 


~Ai / 2 (n, A i) 


(B.25) 


- 66 - 


















Figure 15. Real (blue lines) and imaginary (red lines) parts of the form factors j4 1 / 2 (t, A) (solid 
lines) and A±/ 2 {t,\) (dashed lines), for A = (mz/mh) 2 T ~ 0.52r. The horizontal lines correspond 
to the asymptotic values 2, 4/3 and 0. The vertical lines correspond to r t ~ 0.13, r = 1, r ~ 1.9 
(A = 1), Tb ~ 230 and t t ~ 1300. 


where X\v = = (mz/mh) 2 T\y — 0.52 tw, and analogously A i = m|/(4m-) — 

0.52rj. Note that only the Z/,;/) vector coupling contributes, g] 7 = T^fu )/2 — QjS 2 . The 
form factors are given by 


Ai(r, A) = 2 [3 + 2r - 2A(1 + 2r)] Ii(r, A) - 16(1 - A)J 2 (r, A) , (B.26) 


Ai/ 2 (t, A) = 4 [7 2 (r, A) - Ii(r, A)] 


(B.27) 


where 


h(r,X) 


1 f(r)-f(X) A[g(r)-g(A)] 

2(r-A) 2(r — A) 2 (r - A) 2 ’ 
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The normalisation is chosen to match with the 77 form factors: All (r, 0) = All (r) and 
Ai/ 2 (t,0) = z1i/ 2 (t). The behaviour of ^(rjA) is displayed in Fig. 15, for the relevant 
case A = (m^/m/ l ) 2 r. The VF-boson and f-quark summands in Eq. (B.25) take the value 
Ai(t\v,Xw) — —6.64 and N ct Q t gYA 1 / 2 (T t , Ai)/c^ ~ 0.37, with the lighter SM fermions 
adding a very small contribution. 

In a generic fermionic extension of the SM, we find a decay width 


T(h -> 7Z) 


9 9 *3 / 9 \ 3 
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with the CP-even and odd fermionic amplitudes given by 
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The explicit expression of the four independent form factors will be given below. As far as 
we know, this expression for T(h —> 7 Z), corresponding to a generic set of fermions, was not 
available in the literature. Here the sum runs over all pairs of fermion mass eigenstates: the 
triangular fermion loop is formed by one fj propagator from the h vertex to the Z vertex, 
and two f k propagators from Z to 7 , and from 7 to h. As both h and Z can have off- 
diagonal couplings, j and k can be different. Note that only those combination of couplings 
that are even under the charge conjugation C contribute, because the transition h —> Z 7 
is even: under C, one has Q —> —Q, g^j —*■ —gjj* and x t j —> Xj t = x*-. for x = g A , y , y. The 
P and CP-even (odd) amplitude corresponds to an even (odd) number of axial-vector and 
pseudo-scalar couplings g A and y. 

Let us discuss first the loops involving one fermion mass eigenstate only (j = k). The 
diagonal couplings g^’ A , y% and y, are all real, therefore the form factors bi/ 2 and b\/ 2 are 
irrelevant, while the others reduce to 


ai/ 2 (m, m, m ) = A 1/2 (t, A) , ai/ 2 (m, m, m) = A 1/2 (t, A) = 4J 2 (r, A) . (B.33) 

These two form factors are displayed in Fig. 15 as a function of r. As usual the normalisa¬ 
tion matches with the 77 form factors, in particular A\/ 2 {t. 0) = Ai/ 2 (t). Comparing with 
Eq. (B.10), one finds that the contributions of such fermion ‘diagonal’ loop to the effective 
h'yZ couplings are 


c /i 7 Z 
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Let us now discuss the loops involving two fermion mass eigenstates (j 7 ^ k). The form 
factors are given by 
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Figure 16. Form factors for the triangle loops involving two fermions of mass m and m! = rm , 
for r = 3 and r = 10, in the left and right-hand panels, respectively. We show the real (in blue) and 
imaginary (in red) parts of A 1 / 2 (t, X,r/r 2 , X/r 2 ) (solid lines, upper panels), A 1 / 2 {t,X,t/ r 2 ,X/r 2 ) 
(dashed, upper), S 1 / 2 (r, A,r/r 2 , X/r 2 ) (solid, lower) and B 1 / 2 (t, X,r/r 2 , X/r 2 ) (dashed, lower), as 
a function of r = 4m^/m 2 , for a fixed value of A = (mz/mj) 2 T ~ 0.52 r. The horizontal lines 
correspond to the asymptotic values for the case r = 1 (see Fig. 15). The vertical lines correspond 
to the third family masses, r t ~ 0.13, r& ~ 230 and r T ~ 1300, and to the threshold values, 
m + m' = rrih [r = (1 + r) 2 /4], and m + m' = mz [t = (nih/mz) 2 ( 1 + r) 2 /4]. 
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where Aq, Bq and Cq are the standard Passarino-Veltman scalar functions [122, 123], in the 
convention of Ref. [124], In the literature, these ‘off-diagonal’ loops are often neglected but, 
in models with significant h and Z off-diagonal couplings, they may provide a contribution 
to the decay width of the same order as the ‘diagonal’ loops. For example, the form factors 
in Eqs. (B.35)-(B.37) have been already employed in supersymmetric models to compute 
the charginos loops [125-127]. 

It is useful to combine all the loops involving two given fermions fj and f k . As 
R e (s k jUjk ) = HB(gJ k ykj) and Im (g kj yj k ) = -Im (gf k y k j), the combinations that appear in 
the total amplitude are 


1 


-^1/2 (l~j ) ! Tfc) ^fc) - rj [®l/2 4“ ®l/2 ITlj , Ulj )] 
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= x [h/ 2 (mj,m k ,m k ) - b 1 / 2 (m k ,mj,mj)] , 
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(B.39) 


and analog definitions for the CP-odd counterparts A 1 / 2 and Bi/ 2 . In Fig. 16 we illustrate 
the behaviour of these four form factors as a function of r = jn^/(4m 2 ), where m is the 
mass of the lightest fermion, for a fixed value of the ratio r = m!/m, where ml is the mass 
of the heaviest fermion. In the limit r = 1 one recovers the diagonal form factors, shown 
in Fig. 15. Since w! is the mass of a new charged fermion, it should be sufficiently large 
to comply with experimental lower bounds; requiring for example rrl > m k , one finds that 
only the region r < r 2 /4 is relevant for phenomenology. Note that the behaviour of the 
form factors as r —> 0 is sensitive to the mass ratio: as r increases from 1 to infinity, the 
asymptotic regime settles at larger values of r, and the asymptotic value of the form factors 
Ai /2 and Ai/ 2 tends to zero as 1 /y/r. The form factors Bi/ 2 and Bi/ 2 are zero for r = 1, 
then become of order one as r grows, then tend to zero in the large-r limit. 

B.3 Experimental constraints on the Higgs couplings 

Here we collect the constraints on the Higgs couplings that we use in our analysis. For a 
given Higgs-decay final state a, the LHC measures the signal strength fx a defined as 

_ cr(pp —> h) Br(h —> a) &{‘PP —> h) T(h —> a) 


Pa 


r SM 

1 h 


a SM (pp h ) Br SM (h —> a) a SM (pp —y h) Y SM (h, —> a) 


(B.40) 


where is the total Higgs width. In Table 4 we report the present determination of p a 
for a = 77, ZZ*, WW *, bb, rr, 7 Z, pp as well as on the invisible width, by the ATLAS and 
CMS collaborations. The expected precision for a luminosity of 300 — 3000 fb _1 at 14 TeV 
[128, 129] is reported for these same channels in Table 5. 

From a global fit of the Higgs data, one can also extract information on the Higgs 
coupling to gluons. Taking a rough extrapolation from the fit in Fig. 16 of Ref. [35], we 
find 0.5 < R gg < 1.8 at 99% C.L., where R gg is defined in Eq. (3.4). In the same way we 
extract the analogue quantity for the diphoton channel, 0.5 < i? 77 < 1.9 at 99% C.L., that 
we employ throughout this paper. 
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ATLAS 

CMS 

ATLAS+CMS 


7+8 TeV 

7+8 TeV 

7+8 TeV 

£ 

4.7 + 20.3 fb- 1 

5.1 + 19.7 fb- 1 


m h 

125.36 ± 0.41 GeV [14] 

125.02±g;§? GeV [12] 

125.09 ± 0.24 GeV [13] 

M77 

[14] 

1.12 + 0.24 [12] 

i-i3±8:it [35] 

AAZZ*->4Z 

!-46 toil [14] 

1.00 + 0.29 [12] 

l-29±g;|§ [35] 

M ww* 

l-lS-oJi [14] 

0.83 + 0.21 [12] 

i-os±8:?i [35] 

fJjTT 

1-44 [14] 

0.91 + 0.28 [12] 

1.0718:1 [35] 

A L bb 

0-63±g;P [14] 

0.84 + 0.44 [12] 

0.6518:1 [35] 

M global 

1-18Zuu4 [14] 

1.00 + 0.14 [12] 

1.0918:1 [35] 

W f z 

< 11.0 [130] 

< 9.5 [131] 


r invisible /F h 

< 0.29 [132] 

< 0.58 [133] 



Table 4. Higgs signal strengths measured by the ATLAS and CMS collaborations at y/s = 7 and 
8 TeV and their combination. The error bars and upper limits correspond respectively to ±1 <t and 
95% C.L.. In this paper we adopt the value of the Higgs mass and of the signal strengths from the 
combined fit of the ATLAS and CMS data, reported in the last column. 


The new fermions may or may not affect each of the three factors on the right-hand 
side of Eq. (B.40). Let us discuss first the Higgs production cross section. The dominant 
gluon fusion channel can be modified at leading order only by coloured fermions. The 
weak-vector-boson fusion and associated production can be modified, at tree level, only 
by fermions mixing with the initial state quarks: as we limit our analysis to mixing with 
the third family, we neglect these modifications. Finally, the tt associated production can 
be corrected by those fermions mixing with the top quark. Concerning the total Higgs 
decay width T/,, the dominant branching ratio into bb is affected by fermions mixing with 
the bottom quark, and the second dominant decay channel into WW* is not modified by 
new fermions at leading order. Finally, T^ may be modified significantly by new invisible 
decays, that are possible in the presence of sterile neutrinos. When both a(pp —> h ) and 
T/j are close to their SM value, the signal strength in Eq. (B.40) reduces to p Q ~ R a , where 


R = nh 

±Lry - 


a 


) \^SM + “^neJ 2 + \^n 


T SM {h^a) 


M: 


(B.41) 


SM I 


Here A® ew and are the parity-even and odd new physics amplitudes, respectively. The 
approximation p a ~ R a holds for all colourless new fermions, with the possible exception 
of light sterile neutrinos. 

Several groups performed global fits of the Higgs couplings to the SM particles, allowing 
for deviations in both the fermionic and bosonic decay channels, see e.g. Refs. [134-141]. 
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ATLAS 

[128] 

CMS 

[129] 

a / l^a 

300 fb- 1 

3000 fb- 1 

300 fb- 1 

3000 fb- 1 

77 

0.13 

0.09 

[0.06,0.12] 

[0.04,0.08] 

ZZ 

0.11 

0.09 

[0.07,0.11] 

[0.04,0.07] 

ww 

0.13 

0.11 

[0.06,0.11] 

[0.04,0.07] 

TT 

0.21 

0.19 

[0.08,0.14] 

[0.05,0.08] 

bb 

0.26 

0.14 

[0.11,0.14] 

[0.05,0.07] 

7 £ 

0.46 

0.30 

[0.62,0.62] 

[0.20,0.24] 

PR 

0.39 

0.16 

[0.40,0.42] 

[0.20,0.24] 

r invisible/^ h 

< 0.22 

< 0.14 

< [0.17,0.28] 

< [0.06,0.17] 


Table 5. Expected relative uncertainty at ler on the signal strengths for ATLAS and CMS. The 
expected precisions correspond to yfs = 14 TeV and C = 300 and 3000 fb -1 . We also display the 
expected limit at 95% C.L. on Pinvisible/^ h for the same luminosities. For CMS the two numbers 
correspond to two different estimations of the future uncertainties [129]. 

The fit of Ref. [142] analyzed deviations in the Higgs couplings in the presence of new 
fermions only. 
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